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Abstract 

The  research  contained  in  this  thesis  continues  to  explore  the  c<nicq>ts  Auto¬ 
mated  Formation  Flight  Control  documented  in  three  previous  AFIT  theses.  The  generic 
formation  analysed  consists  of  a  Leader  and  Wingman,  with  the  Wingman  referencing  its 
maneuvers  off  of  Leader  maneuvers.  Specifically,  planar  formation  fiight  control  ccmceming 
only  heading  and  velocity  changes  is  considered.  Next,  the  vertical  separation  constraint  is 
relaxed  to  allow  wing  maneuvers  outside  of  the  flight  plane  of  the  lead  in  order  to  minimise 
the  energy  expended  by  the  wing  in  a  maneuver.  Analysis  of  the  two  forms  oi  formation 
flight  control  investigated  in  this  thesis  reveals  the  close  relationship  between  formation 
geometry,  aircraft  time  constants,  controller  gains,  formation  performance,  and  control 
system  stability.  Integral  control  action  is  determined  to  be  necessary  for  formation  flight 
control.  Nonlinear  simulations  are  accomplished  on  a  digital  computer  to  validate  the  anal¬ 
ysis  of  the  automated  formation  flight  control  syston.  Comparisons  are  made  between  the 
two  forms  of  formation  flight  control  considered,  and  a  third,  energy  conserving  maneuvers, 
in  order  to  determine  which  is  best  for  each  phase  of  flight. 
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AUTOMATION  OF  FORMATION  FLIGHT  CONTROL 


I.  Introduction  To  Formation  Flight  Control  Thesis 
1.1  Overview  of  Thesis 

This  thesis  effort  is  different  than  other  theses  in  the  area  of  flight  control.  Its  purpose 
is  not  to  synthesise  a  controller  to  provide  desirable  flying  qualities  of  an  aircraft.  Instead, 
the  emphasis  is  on  controlling  a  formation  of  aircraft  equipped  with  flight  control  systems. 
This  is  a  high  level  control  problem. 

The  first  chapter  provides  an  introduction  to  the  Automatic  Formation  FUght  Control 
problem  analysed  in  this  research.  It  gives  the  background,  the  assumptions,  the  research 
questioiu  answered,  and  the  scope  of  the  thesis. 

The  second  chapter  reviews  the  literature  pertaining  to  the  formation  flight  control 
problem.  It  gives  a  brief  overview  of  related  work  and  discusses  the  relevance  of  this  thesis. 

The  third  chapter  provides  a  “blueprint*  of  the  model  and  simulation  development 
used  in  this  thesis.  It  is  hoped,  that  on  reading  this  chapter  and  understanding  the 
intricacies  of  the  simulation  models,  one  will  be  able  to  reproduce  the  simulations  performed 
in  this  thesis. 

The  fourth  chapter  begins  the  new  research  into  the  area  of  Automatic  Formation 
Flight  Control  by  extending  Busogany’s  research  from  December  of  1992.  It  is  the  develop¬ 
ment  of  the  analytical  formulation  for  the  formation  flight  controller,  using  two  dimensional 
maneuvers  (i.e.,  the  wing  aircraft  follows  lead  aircraft  maneuvers,  maintaining  the  same 
altitude  as  the  lead).  The  difference  between  last  year’s  approach  and  this  year’s  is  the 
inclusion  of  additional  measurements,  i.e.,  the  heading  error  and  velocity  error,  and  the 
use  of  a  linear  mixer.  This  mandates  the  incorporation  of  the  lead  aircraft  states  into  the 
system  model.  Previously,  the  lead  aircraft  states  were  assumed  to  be  the  disturbances 
affecting/driving  the  formation  flight  control  system.  Now,  the  lead  aircraft  commands 
play  the  role  of  disturbances  driving  the  formation  flight  control  system. 


The  fifth  chapter  extends  the  energy  conserving  maneuvers  concept  first  investigated 
by  Busogany  and  explores  energy  excursion  minimizing  maneuvers.  In  Buzogany’s  thesis, 
the  wing  aircraft  was  allowed  to  moderate  its  speed  by  making  excursions  outside  the  flight 
plane  of  the  lead  aircraft,  in  order  that  the  specific  energy  of  the  wing  be  maintained  at 
a  constant  level.  However,  wing  aircrsit  capability  is  decreased  because  only  formation 
heading  changes  could  be  made.  Now,  the  wing  aircraft  is  still  allowed  to  deviate  from  the 
lead’s  flight  plane,  however,  it  now  can  make  velocity  and  altitude  changes  in  addition  to 
heading  changes,  while  at  the  same  time  minimizing  the  specific  energy  excursions. 

The  sixth  chapter  compares  the  three  concepts  of  formation  flight  control  investigated 
in  this  and  previous  theses.  The  relative  advantages  of  each  method  are  ascertained. 

The  final  chapter  sununarizes  the  finding  of  this  research,  provides  conclusions,  and 
makes  recommendations  for  further  research. 

1.2  Formation  Flight  Contnd  Motivation 

Air  Force  Airborne  Special  Operations  Forces  must  be  able  to  operate  in  hostile 
environments  in  order  to  successfully  complete  their  mission.  In  addition  to  weather, 
terrain  avoidance,  and  enemy  opposition,  their  missions  may  require  formation  flights, 
composed  of  dissimilar  aircraft,  flying  at  low  altitudes  in  close  proximity  to  each  other. 
The  increased  additional  pilot  workload  imposed  on  pilots  in  these  situations  is  known  to 
result  in  pilot  fatigue  and  decreased  mission  effectiveness.  Some  reliable  and  safe  form  of 
automatic  formation  control,  similar  to  an  automatic  pilot,  would  be  beneficial.  Hence, 
the  need  to  decrease  pilot  work  load  is  the  impetus  for  research  into  the  are',  of  Automatic 
Formation  Flight  Control.  This  technology  will  increase  aircraft  safety  during  prolonged 
missions. 

1.3  Background  Information  For  Thesis  Effort 

C-130  aircraft,  models  A  and  B,  are  the  aircraft  models  employed  in  this  research. 
They  are  chosen  for  the  following  reasons: 

1.  C-130s  currently  participate  in  Special  Operations  Missions 
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2.  C-130s  were  uied  in  the  previoiu  three  thetis  efforts  in  this  area,  thus  providing  a 
baseline  by  which  new  research  can  be  compared 


Two  generic  formations  were  investigated  in  previous  thesis  efforts  and  are  reviewed 
in  this  chapter.  These  formations  are  the  trail  and  diamond  formations  illustrated  in 
Figures  1.1  and  1.2,  respectively. 


Figure  1.1  I^ail  Formation  [2:Figure  1.1] 


The  formations  chosen  have  good  operational  characteristics.  The  trail  formation 
reduces  the  probability  of  detection  by  hostile  ground  forces  since  a  minimum  amount  of 
land  is  traversed  by  the  aircraft  in  formation.  In  contrast,  the  diamond  formation  wiU 
over-fly  a  larger  patch  of  the  earth,  however,  it  provides  good  visibility  between  various 
aircraft  [2:page  1-2]. 

A  formation  is  a  not  a  static  entity.  Throughout  the  course  of  a  mission,  a  typical 
aircraft  formation  will  make  several  maneuvos,  initiated  by  the  lead  aircraft,  that  the 
wing  aircraft  must  also  make  while  still  maintaining  the  formation.  Typical  maneuvers 
include  changes  in  velocity,  heading,  altitude,  or  a  combination  thereof.  An  example  of 
these  maneuvers  is  illustrated  in  Figures  1.3  and  1.4,  where  a  diamond  formation  heading 
change  and  a  trail  formation  altitude  change  are  respectively  shown.  In  addition,  the 
formation  geometry  may  not  remain  constant  throughout  the  flight.  It  is  quite  possible 
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that  the  fonnation  geometry  will  change  as  the  threat  changes.  An  example  of  this  is 
shown  in  Figure  1.5  [2:page  1-3]. 


Figure  1.3  Diamond  Formation  Heading  Change  Maneuver 


This  research  is  particularly  concerned  with  the  diamond  formation.  It  was  used 
extensively  by  Buzogany.  Since  Buzogany’s  thesis  is  used  as  a  baseline,  the  diamond 
formation  is  simulated  exclusively,  due  to  time  constraints. 
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Figure  1.5  IVail  To  Diamond  Formation  Change  Maneuver  [2:Figure  1.5] 

1.4  Summary  of  the  Current  Literature 
See  Chapter  n. 

1.5  Research  Objective  and  Questions 

The  thrust  of  the  current  research  into  Formation  Flight  Control  is  to  continue  the 
work  initiated  by  Rohs,  Dargan,  and  Buzogany.  This  entails  enhancing  the  modelling  com¬ 
plexity  of  formation  flight  control,  performii^  the  required  analysis,  and  concept  validation 
that  includes  nonlinear  computer  simulations  of  the  Formation  Flight  Control  System.  The 
foUowing  questions  are  answered  in  this  research  effort: 

1.  In  previous  research,  it  was  assumed  the  lead  aircraft  states,  V'i;  and  1^,  the  heading 
and  velocity  of  the  lead  aircraft,  were  not  measurable  by  the  wing,  and  were  the 
disturbances  driving  the  formation  flight  control  system.  How  is  performance  affected 
when: 

(a)  the  lead  aircraft  states,  and  Vi,  are  observables  incorporated  into  the  for¬ 
mation  flight  control  model 
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(b)  the  lead  aircraft  commands,  V’l.  <uid  are  the  disturbances  of  the  Fbrmation 
Plight  Control  System? 

2.  Can  a  dominant  pole  be  found  through  gain  adjustments  of  the  controller  that  will 
achieve  “good”  formation  response  to  commanded  inputs  (See  Section  1.8  for  defini¬ 
tion  of  “good”)? 

3.  Is  there  a  graphical  way  to  investigate  the  stability  of  two  dimensional  maneuvers? 

4.  Previous  research  investigated  the  concept  of  three  dimensional,  energy  conserving 
maneuvers.  The  concept  is  valid  for  heading  change  maneuvers  only.  It  cannot 
perform  formation  energy  change  maneuvers  adequately.  Is  it  possible  to  enable 
a  formation  to  do  heading,  velocity,  and  altitude  changes  while  minimising  wing 
aircraft  specific  energy  or  excursions? 

5.  Compare  the  formation  flight  control  concepts  of: 

(a)  Two  Dimeiuional  (Planar)  Formation  Flight  Control 

(b)  Energy  Conserving 

(c)  Energy  Excursion  Minimizing 

1.6  Aasumptiona 

The  assumptions  need  to  be  stated  clearly  in  order  to  give  the  reader  a  complete 
picture  of  the  problem.  There  are  many  assumptions  necessary  for  this  problem.  This 
thesis  effort  is  a  building  block.  It  adopts  many  of  the  same  assumptions  used  in  previous 
theses.  In  time,  an  effort  such  as  this  might  result  in  an  operational  system  using  the 
concepts  outlined  in  this  thesis.  If  an  operational  implementation  is  ever  to  be  attained, 
small  stepping  stones  like  this  thesis  need  to  be  placed  in  order  to  provide  a  path  toward 
a  ’’real  world”  system.  The  assumptions  are: 

1.  Each  aircraft  in  the  formation  has  the  following  autopilots  in  place: 

(a)  Mach-Hold  Autopilot 

(b)  Heading-Hold  Autopilot 
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(c)  Altitude-Hold  Autopilot 

All  aircraft  in  the  formation  are  coordinated.  A  coordinated  aircraft  allows  pilots 
to  bank-to-tum  without  the  problems  associated  with  side  slip  velocity  and  adverse 
yaw.  In  other  words,  the  flight  control  system  works  the  rudder  for  the  pilot  in  a 
turning  condition  to  provide  coordination  [3:page  55]. 

2.  No  communication  is  necessary  between  the  aircraft  in  the  formation  in  order  for  the 
formation  flight  control  system  to  work. 

3.  Each  wing  aircraft  possesses  a  perfect  on-board  sensor  capable  of  providing  precise, 
real-time,  position  information  with  respect  to  the  lead  aircraft  in  the  formation.  No 
data  delay  is  assumed. 

4.  Initial  conditions  for  all  simulations  are  straight  and  level  flight  in  a  constant  forma¬ 
tion. 

5.  Wing  aircraft  are  controlled  by  the  automatic  formation  flight  control  system.  They 
automatically  track  maneuvers  made  by  the  lead  aircraft. 

6.  In  order  for  an  analysis  to  be  undertaken,  the  standard  “small  perturbations”  flight 
control  assumptions  must  be  made:  aircraft  can  be  modelled  by  a  linear  system, 
the  earth  is  flat,  etc.  Aircraft  are  highly  nonlinear  systems.  Nonlinearities,  such 
as  kinematics,  induced  nonlinearities  and  saturations  do  play  a  prominent  roll  in 
formation  flight  control.  A  time  invariant  aircraft’s  only  changing  parameters  are 
its  inputs  of  thrust,  and  control  surface  deflection.  Now  the  aircraft  model  will  only 
respond  to  control  inputs  and  not  changing  aircrzdt  parameters.  These  are  valid 
assumptions  because  the  analysis  time  of  interest  is  typically  anywhere  from  1  -  100 
seconds  after  a  maneuver  has  been  performed.  During  this  time  the  center  of  gravity 
is  not  moving  a  noticeable  amount,  nor  does  the  urcraft’s  mass  significantly  change. 

J.7  Scope 

This  research  will  use  first-order  models  developed  by  Rohs  [10]  and  Dargan  [4].  The 
second-order  models  developed  by  Buzogany  wiU  not  be  used  but  will  be  recommended  as 
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an  exteniion  to  the  analyiet  pretented  in  this  thetit.  Both  linear  and  nonlinear  tinmlationi 
are  perfonned  in  order  to  validate  and  *%ie  tone”  linear  analym. 

1.8  Figures  of  Merit 

There  are  no  established  “flying  qualities”  requirements  for  formation  flight  control. 
Safe  maneuvering  in  a  formation  entails  avoidance  of  other  aircraft  or  the  ground.  Beyond 
this,  there  is  a  need  for  wing  aircraft  to  be  able  to  follow  the  lead.  Additionally,  large 
maneuver  command  inputs  can  produce  unstable  responses  due  to  the  nonlinear  saturations 
in  the  system  [2:page  2-14].  These  common  sense  requirements  can  now  be  put  into  more 
technical  specifications: 

•  Automatic  Formation  Flight  Control  System  must  be  able  to  track  commanded  in¬ 
puts  with  sero  steady  state  error  (i.e.,  wing  idrcraft  follow  the  lead  aircraft) 

•  The  system  must  have  “good”  transioit  behavior  (i.e.,  avoid  collisions  between  air¬ 
craft) 

•  The  system  must  be  able  to  handle  large  command  inputs  [4:page  3-36]  [2:page  2-14] 

1.9  Materials  and  Equipment 

A  Sun  Sparc  Station  2  is  the  only  hardware  needed  for  this  effort.  In  addition, 
a  control  system  analysis  software  packi^e  is  needed.  Matlab  4.0  by  Mathworks  is  the 
recommended  software.  This  equipment  is  provided  by  the  Department  of  Electrical  and 
Computer  Engineering  in  the  Navigation  and  Flight  Control  Lab,  Room  131,  and  in  Room 
2001,  Building  641,  Air  Force  Institute  of  Technology.  In  addition,  the  thesis  document  is 
written  in  on  a  Unix  platform  (i.e,  the  Sun  Sparc  Station  2). 

1.10  Conclusion 

This  thesis  effort  helps  continue  the  development  of  Formation  Flight  Control.  By 
alleviating  pilot  workload.  Formation  Flight  Control  will  make  aircraft  missions  more  ef¬ 
fective  and  safer  for  both  aircrew  and  aircraft  participating  in  war  time  missions. 
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II.  Literatun  And  Concept  Review 


This  cluq;>ter  explores  the  literature  pertaining  to  the  concept  of  Automatic  Forma¬ 
tion  Flight  Control.  Its  purpose  is  to  inform  the  reader  on  research  and  concepts  that 
directly  apply  to  this  research. 

B.J  Equation  of  Coriolu 

Similar  to  pilots  flying  in  formation  off  the  lead  aircraft,  an  Automated  Formation 
Flight  Control  system  needs  a  frame  of  reference.  Therefore,  it  is  logical  that  all  relevant 
commands  to  the  formation  be  referenced  to  the  lead  aircraft.  To  accomplish  this  task,  a 
rotating  coordinate  system  is  aflSzed  to  and  centered  upon  the  wing  aircraft. 

With  different  coordinate  frames  moving  relative  to  each  other,  the  Equation  of  Cori¬ 
olis  is  used  to  express  one  frame’s  coordinates  in  terms  r  lother.  Blakelock  states  this 
equation  in  the  following  manner:  "The  motion  of  an  object  as  viewed  from  a  reference 
frame  is  equal  to  the  motion  as  seen  from  the  moving  frame,  plus  the  motion  result¬ 
ing  from  the  relative  angulsur  velocity  of  the  moving  frame  with  respect  to  the  reference 
frame”  [l:page  489].  This  fundamental  mathematical  theorem  is  the  basis  for  any  research 
in  this  area.  It  will  allow  the  Formation  Flight  Control  System  to  resolve  all  the  aircraft 
data  from  the  planes  in  formation  to  the  reference  coordinate  system.  Therefore,  position, 
velocity,  and  acceleration  of  each  aircraft  in  formation  can  be  determined.  Without  this 
information.  Formation  Flight  Control  is  not  possible.  In  equation  form,  the  Equation  of 
Coriolis  is  [l:page  490]: 


=  iZy  +  Wip  X  Rp 


where  (Sc#*  Figure  2.1)  [l:page  490]. 

Ri  the  vector  velocity  of  the  point  in  the  i  reference  frame 

Rp  =  the  vector  velocity  of  the  point  in  question  as  seen  from  the  p  reference  frame 

ufip  =  the  angular  velocity  of  the  p  reference  frame  with  respect  to  the  i  reference 

frame 
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Figure  2.1  Inertial  And  Rotating  Framei  of  Reference  [2:Figare  1.7] 


Rf  is  the  position  vector  of  a  point  in  the  p  frame 

2.2  Iniraformation  Pontioning  System  (IFPS) 

The  Flight  Dynamics  Lab  (WL/FIGS)  is  currently  in  the  process  of  demmutrsding 
a  t3rpe  of  Formation  Flight  Control  System.  This  system  is  a  manual  control  system  with 
the  pilot  of  the  wing  aircraft  flying  a  trajectory  dictated  by  the  flight  computer  [9].  The 
wing’s  trajectory  is  computed  by  determining  how  far  away  the  wing  is  from  the  lead. 
When  this  is  determined,  the  lead’s  trajectory  is  given  to  the  wing  through  the  Heads-Up 
Display  (HUD)  so  the  wing  pilot  can  maintain  formation  by  flying  according  to  informatimi 
displayed  on  the  HUD[9].  There  are  no  provisions  for  automatic  control  of  the  fcwmation 
[9]. 
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i.S  DecoupUug  of  Output:  Porter’*  Deeiqu  Method 

Decoupling  of  output!  in  a  multivariable  control  lystem  is  extremely  desirable.  If  all 
the  states  of  a  system  are  accessible,  output  decoupling  may  be  achieved.  However,  it  is  not 
always  easy  to  gain  access  to  all  of  the  states  of  a  system.  The  feedbsKk  loop  of  any  aircraft 
is  not  a  wire  connected  to  the  output  of  the  aircraft.  On  the  contrary,  sensors  provide  the 
necessary  feedback  in  order  to  close  the  loop.  It  would  be  better  if  output  decoupling 
could  take  place  with  a  smaller  number  of  outputs,  thus  the  inaccessible  states  would 
not  be  needed.  A  control  design  method  which  uses  only  output  feedback  to  generate  an 
error  vector  avoids  the  requirement  for  measuring  or  reconstructing  the  entire  state  vector. 
Dr.  Brian  Porter  developed  a  method  of  designing  a  high  gain  proportional  plus  integral 
(PI)  controller  which  produces  output  decoupling  and  leads  to  very  fast  tracking  of  the 
command  input  by  the  output.  The  standard  MIMO  (multiple-input,  multiple-output) 
plant  is  represented  by  the  following  state  and  output  equations  [2:page  1-12]  [4:page  2-4] 
[S:page  660]: 


x(t) 

=  Ax(t)  -1-  Bu(t) 

(2.1) 

y(*) 

=  Cx(t) 

(2.2) 

where, 

x(f)  =  state  vector 
A  =  plant  matrix 
B  =  input  matrix 
u(f)  =  input  vector 
y(t)  =  output  vector 
C  =  output  matrix 

According  to  the  D’Asso  and  Houpis  classical  textbook  [5:page  660],  this  design  is 
dependent  upon  the  first  Markov  parameter.  The  first  Markov  parameter  is  equal  to  the 
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Figure  2.2  Output  Feedback  IVackiiig  Syitem  With  PI  Ctmtitdltf  [5:Figare  20.1] 


Figure  2.3  Output  Feedback  IVacking  System  With  PI  Controller  Utilising  Measurement 
Matrix  To  Correct  For  Irregular  Plant  [5:Figure  20.4] 


matrix  product  CB.  The  rank  of  this  product  indicates  whether  the  plant  is  one  of  two 
types:  1)  If  CB  has  full  rank,  the  plant  is  considered  regular,  and  2)  otherwise,  the  plant 
is  considered  irregular.  Begular  plants  hare  a  proportional  plus  integral  (PI)  control  law 
implemented  in  the  forward  path  of  the  control  system  operating  on  error  taken  from  the 
input  minus  the  output.  Irregular  plants  have  a  PI  controller  augmented  with  an  inner-loop 
providing  extra  measurements  for  control  purposes  [2:page  1-12]  [4:page  2-4]  [5:page  660]. 
Figures  2.2  and  2.3  show  regular  and  irregular  plants  with  their  respective  implemented 
PI  controllers.  These  Figures  are  taken  from  Figures  20.1  and  20.4  of  D’Asso  and  Houpis 
[5]. 
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t.4  AFITThete* 

Three  AFIT  thetee  have  descrihed  the  topic  of  Formation  Flight  Control.  The  first 
was  presosted  in  March  of  1991.  These  works  are  the  most  current  documented  research 
on  this  topic  [9].  This  thesis  activity  was  sponsored  by  the  Flight  Dynamics  Lab  to  look 
into  the  feasibility  of  Automatic  Formation  Flight  Control. 

Eohs’  thesis  was  the  first  (March  1991).  His  research  lays  the  groundwOTk  for  the 
investigation  of  Automatic  Formation  Flight  Control.  He  condnded  that  manual  flight  con¬ 
trol  system,  such  as  the  IFPS,  could  not  guarantee  formation  maintenance  and  could  result 
in  a  collision  [10:page  100].  Some  form  of  feedback  was  necessary  to  prevent  the  outmaneu- 
vering  of  a  less  citable  aircraft  [lO.'page  100].  Feedback  of  aircraft  errors  with  respect  to 
the  formation  makes  possible  automatic  dissimilar  aircraft  formation  flight. [10:page  100). 

Dargan  (December  1991)  furthered  this  work  by  investigating  different  feedback  pa¬ 
rameters  to  control  the  formation.  He  concluded  that  velocity  and  heading  feedback  could 
be  used  to  keep  the  required  separation  distances  between  aircraft  [4:page  6-3).  The  feed¬ 
back  gains  were  experimentally  determined  and  consisted  of  a  mix  of  separation  and  ma¬ 
neuver  error  [4:page  6-3]. 

Busogany  was  the  last  thesis  student  to  tackle  this  problem  (December  1992).  He 
developed  more  accurate  aircraft  models  and  extended  the  investigation  of  this  problem 
by  looking  at  large  heading  changes  and  how  they  affect  formation  maintenance  [2:page 
8-5].  In  addition,  he  examined  the  case  of  multiple  wing  formations  [2:page  8-5].  He  found 
that  the  Automatic  Formation  Flight  Control  System  performed  well  for  multiple  wing 
aircraft[2:page  8-4].  His  research  is  the  basis  of  the  work  proposed  in  this  document. 

2.5  Conclusion 

The  majority  of  the  work  accomplished  in  Automatic  Formation  Flight  Control  has 
been  research  performed  at  the  Air  Force  Institute  of  Technology.  The  ongoing  Flight 
D]mamics  Lab  effort  concerns  open  loop  control,  while  the  AFIT  research  is  focused  on 
closed  loop  performance.  The  concepts  and  research  outlined  in  this  chapter  (pve  a  brief 
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overriew  the  pertinent  aspects  of  this  research  problem.  'Hie  present  research  will 
c<mtribute  to  this  body  of  knowledge. 
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III.  Model  and  Simidaiion  Development 


This  chapter  provides  the  necessary  information  to  reproduce  the  results  in  this  the¬ 
sis.  It  is  written  to  give  the  reader  a  brief  overview  of  the  simulation,  which  is  implemented 
in  Matlab.  This  chapter  reviews  each  segment  of  the  simulation  and  gives  an  overview  of 
the  actual  implementation.  Control  Strategies  are  also  discussed. 

S.J  Aircmft  Models 

Previous  research  by  Rohs  [10]  and  Dargan  [4]  developed  first-order  aircraft  models 
that  were  incorporated  by  Busogany  [2]  into  his  thesis.  These  models  approximate  the  re¬ 
sponse  of  a  C-130H  aircraft  with  Mach-Hold  ,  Altitude-Hold,  and  Heading-Hold  automatic 
pilots  in  place  [2:page  2-2).  Because  of  the  three  axes  autopilots  working  in  concert,  the 
aircraft  will  generally  behave  like  an  overdamped  first  or  second-order  system  because  of 
the  compensation  each  autopilot  provides.  Three  distinct  channels  are  considered;  velocity, 
heading,  and  altitude.  They  are  given  in  Equations  (3.1)  -  (3.3). 


V{s) 

-L 

Tv 

V.{s) 

V'(-) 

X 

h{s) 

Th 

K(S) 

(3.1) 

(3.2) 

(3.3) 


where. 

Tv  =  time  constant  of  velocity  channel 
=  time  constant  of  heading  channel 
Tft  =  time  constant  of  altitude  channel 

Table  3.1  shows  the  model  time  constants  for  Equations  (3.1)  -  (3.3).  It  also  provides 
the  upper  and  lower  limits  for  the  rate  limiters  employed  in  the  nonlinear  models.  It  ia 
important  to  realize  that  this  table  will  not  exactly  match  Table  2.1  in  Bnzogany’a  thesis. 
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Table  2.1  in  Bitzogany’a  theaia  liata  the  Upper  Altitude  Rate  Saturation  Limit  aa  8.5 
However,  tn  hia  aimulationa  he  implemewted  thia  ae  8  ■^.  Thia  haa  been  determined  through 
inapection  of  hia  aimidation  implementation  tn  Matrix,,  Syatem  Build  and  ita  reaulHng 
plota.  Since  Buzogany’a  theaia  ia  conaidered  a  baaeline  for  thia  theaia  effort,  the  Upper 
Altitude  Rate  Saturation  Limit  ia  implemented  aa  8  in  aimulationa  and  ia  Hated  aa  auch 

in  Table  S.J.  Thia  will  not  match  Table  S.l,  page  S-S  of  Dargan  [4],  nor  Table  2.1,  page 
2-4  of  Buzogany  [2].  A  aimilar  aituation  exiata  for  the  limita  on  the  Velocity  Error  tn  Talde 
S.l.  Theae  were  not  Hated  in  Dargan’a  Table  S.l  or  Buzogany’a  Table  2.1.  The  valuea  are 
taken  from  Buzogany’a  aimulation  of  the  Velocity  modela  (Dargan  only  burned  hia  C-ISOB 
modela  in  hia  theaia  and  the  Upper  and  Lower  Limita  on  the  Velocity  Error  match  thoae  of 
Buzogany).  The  time  responses  of  C-130A  and  B  aircraft  channels  listed  in  Table  3.1  are 
in  Figure  3.2. 


Table  3.1  Aircraft  Rate  Limits  and  Model  Time  Constants  [2:Table  2.1] 


Aircraft 

Parameter 

Lower  Limit 

Upper  Limit 

TV 

TV 

C-130A 

2  sec 

Velocity 

mmm 

msmm 

Heading  Error 

«ec 

8^ 

^  gee 

C-130B 

5  sec 

2  !»ec 

Velocity 

Velocity  Error 

Heading  Error 

Altitude  Error 

TT'M 

ml 

_ _ 

In  addition  to  the  first-order  models  developed  by  Rohs,  Buzogany  developed  second- 
order  models  by  simulating  the  linear  bare  aircraft  state  space  models  and  equipped  these 
with  the  three  basic  autopilots  (Mach-Hold,  Altitude-Hold,  and  Heading-Hold)  [2:page  2- 
4].  The  method  by  which  he  accomplished  this  is  docmnented  in  Appendices  A  and  B  of 
his  thesis  [2].  Because  of  time  constraints  and  to  limit  the  scope  of  this  thesis,  first-order 
models  are  used  in  analysis  and  simulation.  However,  the  inclusion  of  the  second-order 
models  in  the  analysis  and  simulations  is  a  logical  extension  of  the  research  in  this  thesis. 
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C-130A(soUdline: 


time  (sec) 

Figure  3.2  Simulation  of  C-130. 


3.2  Aircraft  Sensor  Measurements 

Feedback  is  an  inherent  pso't  of  any  robust  control  system.  Feedback  entails  mea¬ 
surements  of  the  output  which  are  subtracted  from  the  command  signal  and  allow  an  error 
signal  to  be  created.  Driving  this  error  signal  to  sero  is  the  fundamental  part  of  feedback 
control.  The  feedback  for  an  Automatic  Formation  Flight  Control  system  is  provided  by 
sensor(s)  on  board  each  wing  aircraft.  These  sensors  must  be  capable  of  providing  position 
and  velocity  information  of  the  lead  aircraft  with  respect  to  the  wing  [2:page  2-5]  [4:page 
3-8).  An  actual  sensor  will  have  dynamics  and  a  time  delay  associated  with  its  operation. 
This  will  further  complicate  the  Formation  Flight  Control  System,  by  making  it  more  un¬ 
stable.  In  this  analysis,  the  sensors  are  assumed  to  be  ideal  (i.e.,  no  dynamics,  no  time 
delay,  unit  gain).  Another  problem  with  aircraft  sensor  measurements  is  the  transmission 
of  each  sensor’s  information  to  the  Formation  Flight  Controller.  There  is  already  a  time 
delay  before  the  information  transmitted.  Thus,  the  time  delay  associated  with  data 
transmission  also  contributes  to  the  instability  of  the  system.  This  research  effort  considers 
ideal  sensors  and  no  time  delay  associated  with  the  data  transmission  is  assumed. 

3.3  Formation  Coordinate  System 

This  thesis  uses  the  same  formation  coordinate  system  as  discussed  in  Section  3.4 
of  Dargan’s  thesis  [4]  and  Section  2.3  of  Buzogany’s  thesis  [2].  It  is  restated  here  for 
completeness. 

This  research  uses  two  coordinate  frames: 

1.  Inertial  coordinate  frame 

2.  Rotating  reference  frame  centered  on  the  wing  aircraft 
[2:page  2-6]  [4:page  3-8] 

The  inertial  frame  is  depicted  in  Figure  3.3.  It  has  latitude,  longitude,  and  altitude  as  its 
axes  and  has  a  stationary  origin  [2:page  2-6]  [4;page  3-9]. 

The  rotating  aircraft  reference  frame  is  attached  to  the  wing  aircraft  and  is  depicted 
in  Figure  3.4.  The  wing  aircraft  provides  the  origin  for  the  frame.  The  x  axis  is  in  the 
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Figure  3.3  Inertial  Reference  F^ame  and  Separation  Distances  [2:Figure  2.3] 

direction  of  the  aircraft’s  flight  direction  (i.e.,  aligned  with  velocity  vector  ),  the  y  axis 
points  out  the  right  wing,  and  the  s  axis  points  down,  toward  the  earth  [2:page  2-6]  [4:page 
3-9].  Relative  position  between  the  wing  and  lead  aircraft  is  easily  obtained  using  the  wing 
aircraft  rotating  reference  frame.  In  addition,  the  distance  measurements  provided  by  the 
actual  on  board  sensors  would  provide  information  in  a  similar  manner  as  the  rotating 
reference  frame  [4:page  3-9). 

3.4  Kinematic  Equations 

Mathematical  expressions  must  be  derived  to  determine  the  relative  separations  be¬ 
tween  the  lead  and  wing  aircraft.  This  has  already  been  done  by  Dargan  and  repeated  by 
Buzogany  [4:page8  3-12  to  3-21]  [2:pages  2-6  to  2-10].  It  is  repeated  under  the  assumption 
that  the  reader  does  not  have  adequate  knowledge  of  the  kinematic  development  leading 
to  the  X  and  y  separation  equations.  The  separation  equations  yield  information  that  on¬ 
board  sensors  would  provide  in  a  “real  world”  implementation  of  a  formation  flight  control 
system  [4:page  3-12].  Specifically,  they  provide  the  x  and  y  separation  distances  between 
lead  and  wing  with  respect  to  the  wing  aircraft  {2:page  2-6]. 
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Figure  3.4  Wing’s  Rotating  Reference  F^ame  and  Separation  Distances  [2:Figure  2.4] 

Using  the  Equation  of  Coriolis  (See  Chapter  II),  the  velocity  of  the  lead  with  respect 
to  the  wing  has  been  found  by  Dargan  as  [4:Eq^tion  (3.11)] 


yw  _ 


w 

WL 


(3.4) 


lead  velocity  in  wing’s  frame  wing  velocity  in  wing’s  frame 


where, 

=  velocity  of  the  lead  aircraft  with  respect  to  the  wing,  in  the  wing’s  reference 

frame 


=  angular  velocity  of  the  wing  in  the  wing’s  reference  frame 

^WL  =  position  of  lead  aircraft  with  respect  to  wing  in  the  wing’s  reference  frame 

=  position  of  wing  aircraft  in  its  own  reference  frame 

=  inertial  velocity  of  wing  aircraft  in  its  own  reference  frame 

=  inertial  velocity  of  lead  aircraft  in  its  own  reference  frame 

Figure  3.5  shows  the  basis  for  development  of  the  pertinent  kinematic  equations.  In  Figure 
3.5,  the  following  observations  can  be  made  (See  [2:page  2-9]): 
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0 

0 

0 

,  RZl  = 

y**' 

II 

0 

0 

Vw 

v/ 

0 

>  Vt  = 

0 

0 

0 

(3.5) 


where, 

il>B  =  headixig  error 

Vh  =  inertial  velocity  of  the  lead  aircraft  in  its  own  frame  of  reference 
Due  to  =  0,  Equation  (3.4)  can  be  rewritten  as  [4:Equation  (3.12)] 


=  (3.6) 

All  terms  of  the  right  hand  side  of  Equation  (3.6)  are  known  with  the  exception 
of  .  In  order  to  determine  this,  a  Direction  Cosine  Matrix  needs  to  be  determined 
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to  rotate  into  the  wing  rotating  reference  frame,  that  yielding  .  By  inspection  of 
Figure  3.5,  Dargan  found  the  Direction  Cosine  Matrix  to  be  [4:Equation  (3.8)]: 


coai^a  -sin^ir  0 
sinV*a  co»^B  0 
0  0  1 


(3.7) 


Premultiplying  by  yields  the  following  equation  [diEquation  (3.9)]  [2:Equation 


(2.8)], 


yw  _  ~ 


Vc.sinV'a 


0 


(3.8) 


With  determined.  Equation  (3.6)  can  have  its  variables  substituted  for  the  ex¬ 
pressions  in  Equations  (3.5)  and  (3.8)  to  yield  [4:Equation  (3.13)  and  (3.14)]  [2:Equation 
(2.10)1 


Vi,  cos  ^js 

0 

Vw 

0 

0 

^WL  — 

Visia-ifa 

- 

0 

X 

yw 

— 

0 

+ 

0 

X 

0 

0 

0 

0 

(3.9) 


Vl  cos  ipB 

Vw 

yw  _ 
*WL  — 

- 

- 

0 

0 

0 

0 

(3.10) 


Separating  the  components  in  Equation  (3.10)  yields  the  following  equations: 


i""  =  Vl cost's  +  Wy""  -  Vw 

(3.11) 

=  Vfc  sink's  -  W®’’' 

(3.12) 
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0 


(3.13) 


In  this  research  effort,  x  and  y  separation  distances  are  fed  back  for  onnparison 
purposes.  In  the  simulation  the  separations  distances  are  found  by  integrating  Equations 
(3.11)  and  (3.12).  Since  the  only  angular  rate  used  in  this  simulation  is  there  is  no 
component  of  velocity  in  the  x'*'  direction,  thus,  i'*'  =  0  [4:page  3-16].  Future  research 
efforts  may  want  to  explore  this  area. 

5.5  Initial  CondiUona 

In  this  thesis,  all  simulations  are  performed  with  the  following  initial  conditions: 

•  0°  magnetic  heading  on  inertial  axes 

•  375  ^  formation  flif^t  velocity 

•  500  ft  separation  in  and  directions 

•  500  ft  altitude  (for  three  dimensional  maneuvers) 

These  initial  conditions  are  chosen  because  Busogany  used  these  in  his  simulations  and  his 
thesis  is  the  baseline  by  which  this  research  effort  is  proceeding. 

5.6  Simidation  Overview 

Previotis  research  into  this  topic  at  the  Air  Force  Institute  of  Technology  (AFIT)  has 
been  conducted  using  Matrix,.  Because  of  availability,  convenience,  and  quality,  Matlab 
and  its  program  for  simulating  nonlinear  dynamic  systems,  Simulink,  is  chosen  as  the 
software  for  simulations  performed  in  this  thesis.  This  necessitates  the  need  to  reproduce 
Buzogany’s  results,  where  appropriate,  before  extending  his  research.  A  laige  amount  of 
time  has  been  spent  in  providing  a  quality  simulation  that  not  only  matches  past  results, 
but  is  implemented  in  a  control  analysis  program  that  is  now  the  AFIT  standard. 

The  simulation  is  implemented  as  done  previously  in  Dargan  and  Buzogany.  A  two 
tiered  control  strategy  is  employed.  The  upper  tier  is  for  controlling  the  entire  formation. 
Its  commands  are  issued  by  the  pilot  in  the  lead  aircraft.  The  lower  tier  is  for  controlling  the 
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Figure  3.6  Upper  and  Lower  Tier  Implementation  Strategy 

individual  aircraft  within  the  formation.  The  individual  wing  aircraft  command  originates 
in  the  PI  controller  of  each  wing  aircraft  [4:page  3-11]  [2:page  2-10).  A  pictorial  illustration 
of  this  concept  is  shown  in  Figure  3.6 

Two  simulations  are  developed  in  this  thesis.  The  first  is  a  Two  Dimensional  ma¬ 
neuvers  simulation  in  which  the  wing  aircraft  follows  the  lead  maneuvers  without  leaving 
the  lead’s  flight  plane.  The  second  is  a  Three  Dimensional  maneuvers  simulation  in  which 
the  wing  aircraft  is  allowed  to  leave  the  lead’s  flight  plane  in  order  that  the  excursions  of 
the  wing’s  specific  energy  are  minimised. 
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J.7  Twm>  Dimentumtl  Cimtrolter  Dtwdopmeni 

Previoof  resesrdi  explored  both  PI  control  and  PI  control  oihnaced  with  «  linear 
mixer.  Analyiii  waa  accomplithed,  however,  with  juat  the  PI  ccmtroUer  in  place.  In  thia 
research,  the  aaalytia  ia  extended  to  include  a  PI  controller  and  a  linear  mixer  in  place. 
Both  typea  of  contndler  are  aimulated;  however,  the  thniat  of  the  inveatigaticm  pertaina  to 
the  incluaion  of  a  linear  mixer  in  the  aimulation. 

In  Buaogany’a  analyaia,  the  variables  available  fm  feedback  woe  sq>aration  (x  or  y 
depending  on  the  channel).  In  the  present  analysis,  a  mix  of  separation  and  maneuver 
error  (velocity  error  for  the  c  channel  and  heading  error  for  the  p  channel)  ia  available  for 
feedback. 

S.7.1  PI  ContToUer  and  Linear  Mixer.  A  controller  uses  the  error  between  the 
command  and  the  output  and  dnvea  this  to  sero.  The  error  is  a  mixture  of  separation 
error  (x  or  y  separation  errors)  and  channel  maneuver  error  (heading  or  velocity  errors). 
The  controller  is  the  PI  controller,  which  commands  the  wing  aircraft  to  nullify  the  error. 
Obviously,  only  the  velocity  and  heading  channels  are  controlled  using  this  two  dimensional 
formation  flight  control  system  [2:page  2-12]. 

The  control  law  perturbation  equations  for  the  two  dimensional  case  are  developed 
in  Chapter  IV  and  are  presented  here  as  in  Busogany  [2:Equations  (2.20)  and  (2.21)]: 


Vw,  =  {kyVa  +  *.«)  +  f  (*v Va  +  k,x)  dt  (3.14) 

Jo 

=  *»  +  M)  +  ^  /  *  (315) 

Jo 

where  in  Equation  (3.14), 

X  =  xY  -  X  separation  error  [units:  ft] 

Vs  —  Vi,  •  Vw,  velocity  error  signal  [units:  ^] 
kv  =  velocity  error  mixer  gain  [units:  sec] 

k.  =  X  separation  error  mixer  gain  [units:  1  for  this  research  effort  (See  Ch^ter  IV)] 
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=  s-dunocl  proportkmal  gain  [units:  tec] 

=  c-channel  integral  gain  [units:  see*] 
in  EquatHm  (3.15), 

t  =  ,  y  separation  error  [units:  ft] 

^a  =  ^t-  ifw,  heading  error  [units:  ^] 

=  heading  mot  mixer  gain  [units:  ft] 

kf  =  j  separation  error  mixer  gain  [units:  1  for  this  research  effort  (See  Chiqiter  IV)] 
kgp  =  p-channel  proportional  gain  [units:  ^] 
kyi  =  y-channel  integral  gain  [units:  ^] 

Additionally,  a  control  law  commanding  an  altitude  change  was  developed  by  Dargan  and 
presented  by  Buiogany  as  [2:Equation  (2.17)]  [4:Equation  (4.1)] 


kw,  = 


kfpZ  "f"  h, 


i  fzdt 
Jo 


no  mixer;  only  s  error  available 


where, 

z  =  —  z'*',  s  separation  error  [units:  ft] 

k,p  =  z-channel  proportional  gain  [units:  unit-less] 

kgi  =  z-channel  integral  gain  [units:  ^] 


(3.16) 


3.8  Three  Dimensional  Controller  Development 

Previously,  Buiogany  kept  the  specific  energy  of  the  wing  aircraft  constant  by  allow¬ 
ing  the  aircraft  the  ability  to  deviate  outside  of  the  lead  aircraft’s  fiight  plane  in  response 
to  velocity  changes,  thereby,  maintaining  a  constant  specific  energy  level  (i.e.,  trading  ki¬ 
netic  for  potential  energy).  Chapter  V  extends  Busogany’s  research  by  investigating  if 
allowing  specific  energy  changes  less  than  those  made  in  two-dimensional  maneuvers  can 
have  the  benefit  of  improved  energy  consumption  with  better  ciq)abilities.  This  is  done  by 
formulating  an  altitude  and  velocity  command. 
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5.8.1  PI  Control  Low*  for  Three  Dimenoiouol  Cote.  Btuogany  developed  an 
altitude  command  in  resp<nue  to  x  separation  error.  This  allows  for  the  wing  aircraft’s 
speed  modulation  and  at  the  same  time  energy  conserving  maneuvers  to  take  place  (i.e., 
the  specific  oiergy  of  a  wing  aircraft  ronains  at  its  nominal  value  during  a  maneuver). 

In  energy  excursion  minimising  maneuvers,  the  wing  aircraft’s  oiergy  perturbations 
are  minimised  throughout  the  course  of  a  maneuver  and  are  preferably  ^proximately  sero. 
This  necessitates  the  use  of  a  velocity  control  law  in  addition  to  the  altitude  control  law  to 
minimise  the  swings  in  the  specific  energy  state  of  the  wing  aircraft  and  at  the  same  time 
to  afford  formation  speed  increases.  The  two  control  laws  for  energy  mifiiTnising  maneuvers 
are: 


(3.17) 
dt  (3.18) 

where,  in  Equation  (3.17), 

hxtp  =  ez-channel  proportional  gain  [units:  unit-less] 
hxMi  =  ez-channel  integral  gain  [units: 
and  in  Equation  (3.18) 

ew  =  specific  energy  of  an  2drcraft  [units: 

hw  =  perturbation  away  from  nominal  for  wing  altitude  [units:  ft] 

Tv„  —  time  constant  of  wing  velocity  channel  [units:  sec] 

=  time  constant  of  wing  altitude  channel  [units:  sec] 
g  =  acceleration  of  gravity  [units:  ^^] 

Vo  =  nominal  velocity  of  the  formation  [units:  ^] 

hxtpv  =  zz-channel,  altitude  error  proportional  gain  [units:  ^] 

hxxiv  =  zz-channel,  altitude  error,  integral  gain  [units:  ^;^] 


hw.  =  k*,-*  +  k 


J*xdt 

Jo 


IW.  =  1  ^  hxxpyK  kxtiy  j  K 

\  '  o  L^fcur  J  ^hw  ' 
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h,  =  hi,  ~  kw)  altitude  error  [units:  ft] 


Additionally,  a  velocity  control  law  is  formulated  without  ew  •  This  control  law  is 


vw.= 


fh,dt 

Jo 


(3.19) 


This  development  is  discussed  in  more  detail  in  Chapter  V. 
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IV.  Horizontal  Formation  FHgkt  Control  -  AncJyais 


This  chapter  develops  the  derivation  of  a  linear  and  parameterised  MEMO  plant.  It  is 
used  to  analytically  determine  the  PI  control  and  linear  mixer  gains  for  the  formation  flight 
control  problem.  Analysis  is  constrained  to  the  flight  plane  of  the  aircraft  (i.e.,  velocity 
and  heading  channels).  A  three  dimensional  analysis  is  considered  in  a  later  cluster. 

4.1  Theoretical  Development 

Following  Busogany’s  research,  the  investigation  of  the  linearised  formation  flight 
control  model  proceeds  as  follows: 

1.  Formation  geometry  is  outlined  and  geometric  parameters  are  identified 

2.  Linearized  equations  for  the  aircraft  models  and  aircraft  kinematics  are  developed 
and  nondimensionalized 

3.  The  X  and  y-channels  are  decoupled  and  augmented  using  a  PI  control  law. 

4.  The  augmented  state  equations  are  analyzed  for  steady-state  error  and  stability 

5.  The  stability  envelopes  for  each  channel  parameters  are  plotted  and  analyzed  to  show 
where  acceptable  gains  may  be  picked 

6.  TVial  and  error  process  of  picking  gzuns  exercising  the  nonlinear  simulation 

7.  Varying  gains  and  analyzing  the  poles  and  residues  of  the  pertinent  responses  of  the 
two  channels 

[2:page  5-1] 

Figure  4.1  illustrates  the  formation  geometry  for  a  two  aircraft  formation.  By  in¬ 
spection  of  Figure  4.1,  the  formation  parameters  are  /,  the  separation  distance  {y/xl  -I-  Po), 
and  a,  the  formation  angle  [2:page8  5-1  and  5-2]. 

The  Cartesian  separation  distances  can  be  found  in  Buzogany  [2][Eqs  (5.1)  and  (5.2)] 


Xo  =  I  cos  a 


(4.1) 
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X 


y 


Figure  4.1  Formation  Geometry 


y„  =  isina 


(4.2) 


where, 

Xo  =  nominal  x  separation 

y„  =  nominal  y  separation 

a  =  separation  angle  =  arctaa  j  — »  without  loss  of  generality,  0  <  a  <  | 
Equations  (4.1)  and  (4.2)  relate  the  key  formation  parameters  to  the  separation  distances. 

First-order  aircraft  models  are  used  to  simplify  the  analysis.  In  addition,  the  analysis 
results  can  be  compared  to  results  obtained  by  Busogany  [2].  Equations  (4.3)  -  (4.7)  are 
the  first-order  C-130  aircraft  models  used  in  prior  research  by  Rohs  [10],  Dargan  [4],  and 
Busogany  [2]. 


Vw  = 

- Fur  - Vw. 

(4.3) 

^  = 

1  ,  1  , 

^ 

(4.4) 

4-2 


(4.5) 


II 

II 

1  .  1  . 
- Vt  +  — 

; 

_ —kw  +  -i-hitr 

A|ir  — 

- nw  T  »»Wc 

(4.6) 

(4.7) 


where 

=  the  velocity  time  coiutant  of  the  wing 
=  the  heading  time  constant  of  the  wing 
=  the  velocity  time  constant  of  the  lead 
=  the  heading  time  constant  of  the  lead 

The  nonlinear  kinematic  equations  relating  the  x  and  y  separation  between  the  lead  and 
wing  aircraft  are  given  by  Busogany  [2:Equations  (5.6)  and  (5.7)] 


X  =  Vl  cost's  +  -  Vw  (4.8) 

y  =  V'xsinV’n-t^*  (4.9) 

Because  Equations  (4.8)  and  (4.9)  are  nonlinear,  they  are  linearised  about  a  steady-state 
operating  point  using  small  angle  approximations  and  the  method  of  small  perturbations. 
This  results  in  the  linear  equations  found  by  Busogany  [2;Equations  (5.8)  and  (5.9)] 


X  =  VL  +  i>wyo-Vw  (4.10) 

y  =  Vl^b  -  i>w^o  (4.11) 

Substituting  Equations  (4.1)  -  (4.6)  into  Equations  (4.10)  and  (4.11)  yields  the  equa¬ 
tions  Busogany  fotmd  in  his  research  [2:Equations  (5.10)  and  (5.11)] 
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X 

y 


Vi  +  I — —il>w  +  — /»ina  -  Vw 

Vi.i^L  -  V»w)  - 


1  ,  1,1, 

- i>w  + - Icosa 


(4.12) 

(4.13) 


Simplifying  the  analysis,  Equations  (4.3)  -  (4.6)  and  Equations  (4.12)  -  (4.13)  are  nondi- 
mensionalised  using  the  following  definitions: 


'  ^  +  (<  i<) 

‘  -  V. 


where, 

I  ~  characteristic  length  of  the  nominal  formation  separation 

i  =  char<tcteri8tic  time 

Vo  =  nominal  formation  velocity 

(Buzogany  [2]) 

/,  t,  and  Vo  are  used  to  nondimensionalize  Equations  (4.3)  -  (4.6)  and  (4.12)  -  (4.13). 
Each  variable  in  these  equations  is  nondimensionalized  and  then  substituted  back  into  the 
original  equation  in  order  to  yield  the  final  nondimensional  equation.  The  nondimensional 
variables  are  denoted  by  the  caret  superscript,  “*■  ”.  They  are  as  follows: 


i  =  f 

=  cos  a 

io  =  ^  =  sina 

II 

VfT, 

v„ 

^ 

II 

Yu 

V, 

TV.  =  ^ 

11 

where, 

X  =  nondimensional  x  separation  error 


(4.16) 
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y  =  nondiroeniional  y  tqparation  oror 
Zo  =  nondimensional  nominal  z  separation 
Po  =  nondimensional  nominal  y  separation 
Vw  =  nondimensional  wing  velocity 
Vw,  =  nondimensional  commanded  wing  velocity 
Vi,  =  nondimensional  lead  velocity 
Vi,^  =  nondimensional  commanded  lead  velocity 
=  nondimensional  wing  velocity  time  constant 
=  nondimensional  wing  heading  time  constant 
fy^  =  nondimensional  lead  velocity  time  constant 
=  nondimensional  lead  heading  time  constant 

The  resulting  nondimensionalised  equations  are  as  follows  (4.16).  (The  notation  is 
no  longer  used.  It  u  considered  extrsmeous  notation.  Therefore,  all  remaining 
variables  in  this  section  are  considered  to  be  nondimensional): 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

Equations  (4.17)  -  (4.22)  are  put  into  state  variable  form.  The  MEMO  linearized 
and  nondimensionalized  state  variable  representation  of  the  plant,  inputs,  outputs,  and 
disturbances  are  represented  by  the  following  equations: 


Vw 

V 


sma  ,  sma  , 

*  =  -Vw ifw  + - +  Vl 

1  wr  1  Tr 

- Vw  H - Vw, 

'fVw  '’Vw 

[cos a  ,  COSO  , 

- 1 1  iffw - Vw,  +  i>L 

'’’■^w  j  ’’'^w 

- ’I’w  H - ’I’w, 

Vt  =  -^Vi,-\-^Vi., 

TVl  TVl 

;  1,1, 

Vl  = - ’IfL  +  - 
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X  =  AX^^BU^^TD 

(4.23) 

Y  =  CX 

(4.24) 

therdbre, 


X  = 


_  _ 

■  0 

-1 

0 

i 

Vw 

0 

_ 1_ 

0 

y 

0 

0 

0 

0 

0 

0 

Vj, 

0 

0 

0 

.  . 

0 

0 

0 

Y  = 


»ma 

^■*w 

0 


Q  una 
1 

0 


0 

0 

0 


1 

^*w 

0 


,  nag 
’••*r 
0 


te-') 


1 

’•♦ir 

0 


0 

0 


0 

0 


1  0  ■ 

c 

0  0 

Vw 

0  1 

y 

0  0 

i>w 

0 

Vl 

0 

^*L  -1 

.  . 

Vw.  ' 

0  0 

■  Vt.  ■ 

+ 

0  0 
_L.  n 

.  . 

0 


X 

■  1 

0 

0 

0 

0 

0  ■ 

0 

-1 

0 

0 

1 

0 

y 

0 

0 

1 

0 

0 

0 

.  . 

0 

0 

0 

-1 

0 

1 

1 

z 

Vw 

y 

i>w 

Vl 

L 


(4.25) 


(4.26) 


where, 


X  =  state  vector 


A  =  plant 

r  =  disturbance  input  matrix 
D  =  disturbance  input  vector 
Y  =  output  vector 
C  =  output  matrix 


4-6 


There  is  almost  total  decoupling  between  the  x  and  y-channel.  The  only  coupling 
between  the  two  occurs  due  to  the  disturbance  in  the  x-channel.  This  is  shown  in  Sec¬ 
tion  4.1.2. 

Due  to  decoupling,  the  state  vector,  can  be  partitioned  as  follows: 


(4.27) 


4.1.1  Y-Channel .  Equations  (4.25),  (4.26),  and  (4.27)  yield  the  decoupled  y-channel 
state  equation: 


y 

^>w 

= 

.  . 

0  (^-i) 

» 

0  0 


1 

0 

I 

■'*L 


y 

•  co»a  ■ 

0 

V’w 

.  . 

+ 

1 

’•♦tr 

0 

+ 

0 

1 

-  ^*1  . 

(4.28) 


The  y-channel  control  variable  is  i^he  wing  aircraft’s  heading  conunand.  It  responds 
to  formation  separation  change  in  the  y-channel.  Note  that  the  y->jiannel  control  system 
now  is  augmented  to  include  the  new  state  variable,  ipi,  which  was  previously  treated  as 
a  disturbance.  Hence,  at  this  level  of  modeling,  an  additional  output  is  now  available  for 
feedback  and  employing  a  mixer  as  in  Dargan’s  thesis  is  needed.  The  role  of  the  disturbance 
signal  is  relegated  to  Mixer  parameters  are  used  to  form  the  generalised  y-channel 
error  signal: 


(4.29) 

where, 

V’e  =  V’l  -  V’w  =  heading  angle  error  (radians) 
ky  =  y-separation  error  signal  mixer  gain 
=  heading  error  signal  mixer  gain 
A  PI  control  law  operating  on  the  error,  e,,  is. 


e,  =  *,y  + 
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where, 


V^.  =  e»  A 


(4.30) 


kgp  =  non-dimeniionai  proportional  gain  operating  on  the  error,  e, 

=  non-dimensional  integral  gain  operating  on  the  error,  e. 

By  differentiating  Equation  (4.30),  the  following  equation  is  obtained, 

(4-31) 

Substituting  Equation  (4.29)  into  Equation  (4.31)  yields, 

(♦•32) 

The  augmented  state  equations  incorporating  both  the  linear  mixer  and  the  PI  controller 
are 


X,  = 

(4.33) 

n  = 

C7,jr, 

(4.34) 

Therefore, 

1  y  ]  f  ° 

ti>w  _  0 

ft  0 

.  kfik. 

fe->)  ■ 

» 

0 

CMO  Hr  1 

-1— 

’•♦tr  , 

0  Vl 

c,,.  J  L  J 

0 

+  1.  V-i.  (4-35) 
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where, 


10  0  0 
0-110 


v 

y 

.  . 

(4.») 


c,,  =  L-t,  -  l)  +  -  V*. 

(4.37) 

+  bmiht 

(4.38) 

«<>•“  +  *#) 

(4.39) 

and 


Xy  =  augmented  y-channel  state  vector 
Af  =  augmented  y-channel  plant 
r,  =  augmented  y-channel  input  matrix 
=  augmented  y-channel  input  vector 
Yf  =  augmented  y-channel  output  vector 
Cf  =  augmented  y-channel  output  matrix 


4. 1.1.1  Static  Stability  Analyaia,  Y-Channel.  It  is  imperative  that  the  wing 
aircraft  be  able  to  track  the  lead  aircraft  with  zero  steady-state  error.  This  requirement 
can  be  stated  as  follows: 


i4s  t  -»  00,  y(oo)  =  0,  and  V^v(oo)  =  ^l(oo)  =  t^.(oo)  =  =  D,  =  1  (4.40) 
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This  reqairement  is  met  by  the  y-rhsamel.  This  is  Uhutreted  by  observing  the  y- 
channel  response  to  a  unit  step  disturbance, 


0 

1 

1 

1 


(4.41) 


The  results  of  Equation  (4.41)  are  true  if  and  only  if  >4,  is  invertible.  This  conditkm 
implies  the  following  relation: 


0  (4.42) 

Conclusion:  hyi  /  0  and  by  ^  0.  Hence,  integral  control  action  is  necessary.  Further¬ 
more,  y-feedback  is  necessary  and  if),  feedback  alone  will  not  suffice.  TherefiDre,  without 
loss  of  generality,  assume  by  =  1.  See  Section  4.2  for  the  rationale  to  allow  by  =  1. 

4. 1.1. 2  Dynamic  StabOUf  Analysts,  Y-Channd.  The  dynamic  stability 
analysis  consists  of  1)  obtaining  the  y-channel  characteristic  equation,  and  2)  performing 
a  Routhian  stability  analysis.  The  characteristic  equation  is 


det(5/  -  i4y)  =  s^  os’  -h  is’  -i- cs  +  d  (4-43) 

where  the  coefficients  are. 


a  = 

(1  +  byy[coBa  -I-  b^])  -I- 

(4.44) 

6  = 

+  “»«])  +  ~ 

(4.45) 

c  = 

^  1  b 

byi - -I- - 

(4.46) 

d  = 

“TSw-rsL 

(4.47) 

The  following  relations  result  from  the  Routh  criterion  application  of  Table  4.1: 
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Tkhle  4.1  Y-Chaimd  Rogthian  Axraj 


4 

a 

a 

1 

0 


1  id 

a  e 

sirs  d 

*m‘d 

d 


KHi* 


row 


;^(l  +  *„cosa  +  *„*^)+:^  >0 


(4.48) 


[(1 +  *»  co«  «  + *wM  (*»  +  *»<*♦  + *r  «)] 

(1  +  k„coia  +  k„k^)*  +  71^  (1  +  k^cosa  +  >  0 


(4.49) 


row  a^ 


(1  +  *„cosa  +  k„k^)  +  (*w  +  +  *r  coaa)  +  >  0  (4.60) 


row  a” 


— —k^  >  0  (4.51) 

Similar  to  Btuogany’s  thesis,  these  stability  conditions  show  the  relationships  be¬ 
tween  controller  and  mixer  gains,  formation  geometry,  and  first-order  model  time  constants. 
In  order  to  guarantee  stability,  all  of  these  inequalities  must  be  satisfied. 
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Aa  in  the  atatic  atability  analyaia,  the  djmamic  atability  analyaia  reveala  the  need 
for  integral  control  action  in  (»d»  to  aatiafy  atability  conditiona.  Remember,  for  theae 
atability  rdationa,  =  1  and  aa  auch  ia  not  aeen  in  the  atability  eqnationa.  The  raticmale 
for  thia  ia  given  in  Section  4.2. 

4.1.2  X-Channel .  Once  again,  Equationa  (4.25),  (4.26),  and  (4.27)  yield  the 
decoupled  x-channel  atate  equation: 


where, 

=  V'w,  -  V'w 

There  are  two  diaturbancea  in  the  x-channel,  the  lead  aircraft  *a  velocity  command, 
Vi,  and  aome  coupling  of  the  y-channel  into  the  c-channel  due  to  heading  changea,  d^. 
Viv,,  the  wing  aircraft ’a  velocity  command,  reaponda  to  formation  aeparation  change  in 
the  x-channel.  Mixer  parametera  are  used  to  form  the  generalised  y-channel  error  aignal: 

e.  =  *,*  +  k,Ve  (4.53) 

where, 

Ve  =  Vi,-Vw  =  velocity  error 
kf  =  x-separation  error  aignal  mixer  gain 
kv  =  velocity  error  aignal  mixer  gain 
A  PI  control  law  operating  on  the  error  aignal,  e.,  ia, 

V  =  +  (4.54) 

Jo 


4-12 


where, 


=  aon-dimeofional,  x-ieparation  error,  proportioui  gem 
=  non-dimemional,  z-aeparation  error,  integral  gain 
By  differentiating  Equation  (4.54),  the  following  equation  is  obtained, 

Vw.  =  +  W  +  k^ik,x  +  k,V,)  (4.56) 

The  augmented  itate  eqtiationa  incorporating  both  the  PI  controller  and  the  linear  mixer 
are 


X.  =  A,X,-\-T.d, 

Y,  =  C,X. 


(4.56) 

(4.57) 


therefore. 


X 

Vw 

Vl 

.  Vw.  . 

0 

0 

0 

kgikg 


-1 
_ i_ 

0 

('Vw 


Y.= 


1 

0 

1 

‘^L 

n 


0 

1 

^w 

0 


X 

Vw 

Vl. 

+ 

Vw. 

0  0  0 
-110 


0 

1 

k  k.— 

**F**rvj 


■iao 

^*w 

0 

0 

L  L  liao 
*’*^’^*W 


X 

Vw 

X 

Vl 

K 

.  . 

Vl. 


(4.58) 

(4.59) 


where. 


(^Vw 

('Vi 


—  kgmkx  “f"  kgpk^  kgikjf 

—  kgmkg  —  kgmk^—  +  kgik„ 


(4.60) 

(4.61) 
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L 


and 


Xg  =  augmented  x-channel  state  vector 
A,  =  augmented  x-channel  plant 
r,  =  augmented  x-channel  input  matrix 
D,  =  augmented  x-channel  input  vector 
Yg  =  augmented  x-channel  output  vector 
Cg  =  augmented  x-channel  output  matrix 


4. 1.2.1  Static  Stability  Analysis,  X-Channel.  Because  it  is  imperative  that 
the  wing  be  able  to  track  the  lead  aircraft  with  sero  steady-state  error,  the  conditions  of 
Equation  (4.40)  can  be  restated,  using  x-channel  parameters  as. 


As  t -*  oo,y(oo)  =  0,  and  Viir(oo)  =  Vi;(oo)  =  Viir,(oo)  =  Vi.  =  1  (4.62) 


The  requirement  of  Equation  (4.62)  is  met  by  the  x-channel.  This  is  illustrated  by 
observing  the  x-channel  response  to  a  unit  step  disturbance, 


*(0„  =  ima-(5/-  >i.)-‘r.  =  -a;^t^  = 


!‘- 

mg 

0 

ling 


(4.63) 


Since  it  is  proven  in  the  y-channel  analysis  that  ^  -*  at  steady-state,  the 
disturbance  due  to  the  coupling  term  -  ^)  will  — »  0  in  steady-state.  Therefore, 

the  second  column  in  equation  (4.63)  does  not  affect  the  response  in  steady  state. 


The  results  of  Equation  (4.63)  are  true  if  and  only  if  A^  is  invertible.  This  condition 
implies  the  following  relation: 
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(4.64) 


Conclution:  ik^  /  0  and  k,  ^  0.  Hence,  integral  control  action  it  necessary.  Further¬ 
more,  c-feedback  is  necessary  and  V’.  feedback  alone  will  not  suffice.  Therefore,  without 
loss  of  generality,  assume  k,  =  1.  See  Section  4.2  for  the  rationale  to  make  k,  =  1. 

4. 1.2.2  Dynamic  Stability  Analysis,  X-Channel.  The  cbaracteristic  equa¬ 
tion  is 

det(5/  -  A.)  =  +  es®  +  fs^-\-g»  +  h  (4.65) 

where  the  coefficients  are, 


- (1  -1-  kv)  +  — 

(4.66) 

-—(Asp  +  *«k,)  -1-  — 

(4.67) 

— +  — 

''Vvr  "^Vv 

(4.68) 

(4.69) 

Table  4.2  X-Cbannel  Routhian  Array 


I* 

.0 


1 

e 


«/-s 

e 


/  h 

9 

h 


The  following  relations  restilt  from  the  Routh  Criterion  application  of  Table  4.2. 
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row 


’Vi 


(4.70) 


row  a’ 


—  [(1  +  k^k,)ik^  +  kM  +  —(1  +  +  j-(l  +  *.,*.)  -  *«  >  0  (4.71) 

’Vw  L  TVt  J 

row 


Tv^  +  TVi(l  +  kgpk^)  +  Ty^(kgf  +  kgik,)  +  Ty^kgi  >  0  (4-72) 


row 


;rV*«  >  ° 

As  in  the  y*channel,  these  stability  conditions  show  the  relationships  between  con¬ 
troller  and  mixer  gains,  formation  geometry,  and  first-order  model  time  constants.  In  order 
to  guarantee  stability,  all  of  these  inequalities  must  be  satisfied. 

As  in  the  static  stability  analysis,  the  dynamic  stability  analysis  reveals  the  need  for 
integral  control  action  in  order  to  satisfy  stability  conditions.  Remember  for  these  stability 
relations  that  jb,  =  1  and  as  such  is  not  seen  in  the  stability  equations.  The  rationale  for 
this  is  given  in  Section  4.2. 

4. 1.2.3  Comparison  of  Analysis  Wiih  Prior  Work.  In  order  to  ensure 
the  previous  analysis  is  correct,  it  is  checked  with  previous  results  obtained  by  Busogany 
[2:pages  5-9  and  5-12].  Busogany’s  Routhian  analysis  did  not  include  a  mixer  in  either 
channel  and  assumed  the  lead  aircraft  states  as  disturbances.  The  analysis  in  Sections 
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4.1.1  and  4.1.2  includes  a  mixer  and  incorporates  the  lead  aircraft  states.  Therefore,  to 
check  the  current  Bx>uthian  analysis  with  prior  analysis,  the  following  conditions  must  be 
applied  to  the  inequalities  of  Ekjuations  (4.48)  -  (4.51)  and  Equations  (4.70)  -  (4.73): 

1.  Turn  the  mixer  off  (i.e.,  set  k,  =  1,  =  0).  Therefore,  e,  =  y 

2.  Lead  aircraft  states  were  not  available  for  analysis  in  the  previous  thesis.  Therefore, 
terms  that  include  the  lead  aircraft  states  are  neglected.  This  is  the  same  as  not 
having  sensors  that  measure  lead  aircraft  parameters. 

After  applying  the  previous  conditions  to  the  pertinent  equations,  the  resulting  in¬ 
equalities  are  the  same  as  those  found  earlier  in  Buxogany  [2:Equations  (5.36)  -  (5.39)  and 
Equations  (5.53)  -  (5.55)].  Therefore,  it  is  shown  that  the  inequalities  derived  in  the  this 
thesis  are  correct  and  match  the  results  of  previous  work. 

4-2  Gain  Envelopes 

The  y-channel  inequalities  in  Equations  (4.48)  -  (4.51)  and  the  x-charnel  inequalities 
in  Equations  (4.70)  -  (4.73)  can  be  employed  to  graphically  portray  the  resulting  stability 
envelope  for  the  respective  gains  of  each  channel.  Tbe^e  are  four  gains  in  each  channel 
(y-channel:  k^,k„,  k^,  and  k,;  z-channel:  k,i,kgp,  k,,  and  k.)  which  influence  stability. 
In  order  to  simplify  analysis,  k,  =  k,  =  1.  This  can  be  done  without  loss  of  generality.  To 
illustrate  this  point,  recall  Equation  (4.29) 

e,  =  k,y  -f-  k^V’e  (4.74) 

kp  can  be  divided  from  both  sides  to  yield, 

^=y-f-^V’«  (4.75) 
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Equation  (4.75)  can  now  be  substituted  into  the  equation  for  the  PI  controller,  Equation 
(4.30) 


""  lo 

The  gains  divided  by  i,  in  the  previous  equations  (Equations  (4.75)  and  (4.76))  can  be 
redefined  to  include  in  their  definitions.  Thus 


(4.77) 


Thus  allowing  ib,  =  1  can  be  done  without  loss  of  generality,  because  the  gain  is  ab¬ 
sorbed  into  the  remaining  gains.  This  same  example  can  be  employed  using  the  s-channel 
equations  to  show  kg  =  1  without  loss  of  generality. 

The  eight  inequalities  in  Equations  (4.48)  -  (4.51)  and  (4.70)  -  (4.73)  are  difficult  to 
plot.  Simplifying  the  inequalities  by  employing  new  variables  via  a  variables  transformation 
will  make  the  inequalities  more  tractable. 

Consider  the  y-channel  first.  The  inequalities  representing  stability  conditions  have 
been  kept  in  factored  form  in  order  that  the  variables  transformation  be  discerned  easily. 
Upon  examination  of  the  inequalities  given  by  Equations  (4.48)  -  (4.51),  it  is  readily  ap¬ 
parent  that  all  four  stability  conditions  have  three  common  factored  expressions.  These 
are  suitable  choices  for  the  proposed  variables  transformation.  Hence,  the  resulting  trans¬ 
formed  variables  for  the  y-channel  aie  defined. 


k  ^ 

1  +  k„  cos  a  +  k„k^ 

(4.78) 

1  ^ 

k„  -1-  k^k^  ■>r  k,i  cos  a 

(4.79) 

k  — 

kfi 

(4.80) 
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A  MTnilaf  examination  it  performed  upon  the  inequalitiet  the  s-channel.  The 
transformed  variables  for  the  c-channel  are: 


m 

1  +  *.,*. 

(4.81) 

n 

A 

-l- 

(4.82) 

A 

b.i 

(4.83) 

The  resulting  inequalities  with  Equations  (4.78)  -  (4.82)  substituted  into  the  original  fac¬ 
tored  inequalities  given  by  Equations  (4.48)  -  (4.51)  and  (4.70)  -  (4.73)  are: 


y  row 


y  row 


1  .  1  « 
- *-»- - >  0 


ki + + 4-*  -  >  0 


(4.84) 


(4.85) 


y  row  S'* 


>  0 


(4.86) 


y  row  s” 


>  0 


(4.87) 


Similarly, 


X  row  <•' 


1  1  o 

- m  H - >  0 

’'Vw  'TVi 


(4.88) 


X  row  s* 


1  1  a  1  .  « 

— -mn  -f - mr  +  -j-m  -  *«  >  0 

‘^Vw'TVi  Ty^ 


(4.89) 
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z  row 

>  0  (4.90) 

y  row  4° 

*^  >  0  (4.91) 

The  traoBformation  of  variables  for  the  y  and  z-channeb  is  a  nonlinear  transformation 
in  which  Ky  =  {h„,  k^)  €  W®,  is  mapped  into  Gy  =  (k,  /,  *^)  G  ^2®,  and  ifjr  = 

{k^,kgi,k,)  €  72^  is  mapped  into  subspace  Gx  =  (nt,n,ftsj)  G  These  nonlinear 
transformations  map  the  respective  gains  into  a  space  in  which  the  stability  equations 
become  tractable.  Once  these  stability  envelopes  are  plotted,  any  choice  of  gains  can  be 
transformed  into  the  new  space  to  determine  if  the  guns  result  in  a  stable  linear  system. 

The  transformed  y  and  z-channel  inequalities  have  a  very  similar  form.  The  sym¬ 
metry  between  the  two  channels’  stability  equations  shows  the  elegance  of  the  modelling 
employed.  The  differences  between  the  two  consist  of  the  pertinent  time  constants  and 
gains  of  the  respective  channel.  This  simplification  allows  easier  parametric  analysis  of  the 
stability  of  the  two  channels. 

The  transformed  inequalities  in  Equations  (4.84),  (4.86),  and  (4.87)  and  their  cor¬ 
responding  transformed  z-channel  inequalities,  are  all  equations  of  planes  in  a  three  di¬ 
mensional  space.  The  inequalities  of  Equations  (4.85)  and  (4.89)  are  equations  for  quadric 
surfaces  [6:pagel31].  The  z  and  y-channel  inequalities  are  plotted  in  Figures  4.2  and  4.4. 
These  figures  show  all  of  the  inequalities  plotted  on  one  plot  with  its  associated  projection 
in  all  three  planes.  Figures  4.3  and  4.5  show  the  stability  envelopes  of  the  y  and  z-channel 
with  all  other  conflicting  constraints  removed.  Notice  that  for  stability,  the  gains  must  be 
between  the  two  surfaces.  Therefore,  it  is  now  possible  to  pick  gains,  parameterize  them, 
and  determine  if  they  yield  a  stable  control  system  by  plotting  them  in  the  respective  gain 
envelope.  The  problem  of  attaining  a  stable  system  has  been  solved.  (Appendix  C  has  the 
individual  plots  of  Equations  (4.84)  -  (4.91)) 
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The  relation!  found  by  Buiogany  and  presented  in  Cluster  5  of  his  thesis,  are  also 
plotted  in  Figures  4.6  and  4.7  (See  Busogany  [2:Equatioxu  (5.36)  -  (5.39)  and  (5.53)  • 
(5.55)]).  The  re^ons  of  stability  in  these  plots  are  shaded  for  ease  of  identification. 


All  Constraints  (xi  Same  I^ot 


K-L  Plane 
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K_yi  -  L  Plane 


K_yi  -  K  Har» 


Figure  4.2  Nondimensional  Y-Channel  Constraints 


Figure  4.3  Nondimensional  Y-Chatmel  Envelope:  Stability  Between  the  Surfaces 
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All  Constraints  on  Same  Plot 


M-NPIane 
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Figure  4.4  Nondimension&l  X-Chsnnel  Constrsints 


The  Sudnlity  Eovelope  for  die  X-Chamel  from  Anodm  Angle 


5  5  N 
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Figure  4.5  Nondimensional  X-Channel  Envelope:  Stability  Between  the  Surfaces 
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Figure  4.6  Nondimensional  Y-Channel,  Buxogany’s  Case:  No  Mixer,  a  =  45”;  Shaded 
Region  Denotes  Stability 
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4.i-l  Bnzogany’s  Teat  Case  Piotted  Innde  SiaUtUy  Envelope.  The  preceding 
figures  show  the  stability  envelopes  for  both  BuK^any’s  case  (PI  controller  only,  Figures 
4.6  and  4.7)  and  the  new  case  (PI  controller  and  linear  mixer,  Figures  4.3  and  4.5).  In  order 
to  insure  stability  envelopes  have  been  properly  defined,  the  test  case  used  in  Busogany’s 
thesis  is  plotted  inside  both  the  two  dimensional  envelope  and  new  three  dimensional 
envelope.  Table  4.3  defines  this  test  case.  Note  that  the  values  listed  are  dimensional. 
However,  the  values  plotted  have  all  been  nondimensionalised. 

Using  the  same  gains  Btuogany  used  in  his  preliminary  research  (see  Table  4.3),  the 
following  relations  for  the  transformed  variables  are  determined  (remember,  no  mixer  was 
used  in  the  previous  analysis,  i.e.,  =  k,  =  0): 


ka  —  1  +  cos  a 
Ib  =  k„  +  k^coaa 

kjfia  —  ^ 


(4.92) 

(4.93) 

(4.94) 


and, 


nifl  =  1 

ng  =  kgf 


k  ■  —  k  ■ 

*X«B  —  *«» 


(4.95) 

(4.96) 

(4.97) 


By  substituting  the  values  for  the  test  case  into  Equations  (4.92)  -  (4.97),  values  can 
be  found  which  are  then  plotted  inside  the  stability  enveloped  developed  earlier.  Plotting 
these  values  yields  the  plots  in  Figures  4.8  -  4.10.  The  values  used  for  the  3-dimensional 
plots  are  shown  in  Table  4.4. 
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Tkble  4.4  Values  of  Test  Case  Plotted  in  S-Dimensional  Nondimensional  Stability 
Envelopes 


PARAMETER 

NONDIMENSIONAL  VALUE 

k 

501 

1 

716.53 

•hi 

13.333 

m 

1 

n 

1.4142 

k,i 

.026667 

Figures  4.8  -  4.10  show  that  using  Busogany’s  gains  yields  a  point  within  the  stability 
envelope,  as  expected.  Variations  are  explored  in  the  next  section. 


kjp  k-V 

Figure  4.8  Buzogany’s  Case  (No  Mixer):  Nondimensional  Stability  Envelopes  With  Test 
Case  Plotted  (See  Figure  4.6  to  Examine  Entire  Buzogany  Y-Channel  Enve¬ 
lope) 


4.2.1. 1  Varying  and  i,.  Buzogany  assumed  k^  =  0.  By  varying  k^,  it  is 
useful  to  see  what  happens  in  the  stability  envelopes,  k^  and  ib,  were  varied  in  a  manner 
according  to  Table  4.5.  The  resulting  plots  are  in  Figures  4.11  and  4.12.  As  can  be  seen  k^ 
and  b,  can  be  varied  and  still  be  within  the  stability  envelope  (The  X-Channel  envelope 
in  this  plot  looks  different  due  to  the  magnitude  of  k^  in  the  previous  plots). 
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SOLI  7liM 


Figure  4.9  Nondimensional  Y  Stability  Envelope  With  Test  Case  Plotted  (Both  Plots  are 
Zoomed  In  To  The  Test  Case  Location:  Right  Plot  Also  Has  K^i  Axis  Scaled. 
See  Figure  4.3  To  See  Whole  Envelope  Plot) 


Figure  4.10  Nondimensional  X  Stability  Envelope  With  Test  Case  Plotted  (Both  Plots 
are  Zoomed  In  To  The  Test  Case  Location:  Bight  Plot  Also  Has  K^i  Axis 
Scaled.  See  Figure  4.5  To  See  Whole  Envelope  Plot) 


.Mixer  Caii 

u  for  Plots 

PARAMETER 

UPPER 

LOWER 

.0001 

-.0001 

.05 

-.05 

By  exploring  the  envelope  in  this  way,  a  designer  is  allowed  to  pick  stable  giuns  and 
check  to  see  where  they  lie  inside  the  envelope. 


ii#‘ 


Figure  4.11  Nondimensional  Y  Stability  Envelope  With  Test  Case  Plotted  and  Varied 
(X  Is  Associated  With  -k^  and  O  Is  Associated  With  k^... Right  Plot  is 
Zoomed  In  Further.  See  Figure  4.3  To  See  Whole  Envelope  Plot)) 


4.2.2  Gain  Envelope  Conclusion.  The  stability  problem  of  the  z  and  y-channels 
is  now  solved.  Picking  gains  to  obtain  a  stable  system  can  be  accomplished  easily.  However, 
stability  of  the  system  is  only  one  aspect  of  the  control  system  design.  Within  these  stability 
envelopes,  there  are  many  combinations  of  gains  that  will  yield  a  stable  system.  However, 
system  performance  in  each  of  these  channels  has  not  been  addressed.  The  performance 
aspect  is  the  true  focus  of  the  research  and  can  now  be  explored  with  the  knowledge  of  the 
stability  boundaries. 


Figure  4.12  Nondimensional  X  Stability  Envelope  With  Test  Case  Plotted  and  k.  Var¬ 
ied  (X  Is  Associated  With  -k,  and  0  Is  Associated  With  ik,... Right  Plot  is 
Zoomed  In  Further.  See  Figure  4.5  To  See  Whole  Envelope  Plot)) 

A  nice  result  of  picking  gains  and  plotting  them  inside  the  stability  envelope  is  an 
indication  of  system  robustness  to  changes  in  the  gains.  A  point  that  is  in  the  middle  of  the 
envelope  can  move  about  the  space  defined  by  the  envelope  and  still  yield  a  stable  system. 
However,  a  point  that  is  near  one  of  the  stability  envelope  edges  is  not  very  tolerant  of 
gain  variations.  A  variation  in  the  wrong  direction  will  result  in  an  unstable  system. 

The  analysis  accomplished  on  gain  envelopes  for  the  Automatic  Formation  Flight 
Control  System  modelled  in  this  chapter  is  very  useful  when  picking  gains  that  yield  good 
performance  emd  robustness.  However,  this  analysis  is  performed  in  the  linear  domain 
only.  It  is  invalid  when  the  system  is  nonlinear. 

4-3  Pole  Placement  Through  Controller  and  Linear  Mixer  Gain  Adjustment 

The  closed-loop  MIMO  state  equations  developed  in  Sections  4.1.1  and  4.1.2  are  used 
to  aid  in  choosing  suitable  gains  for  the  PI  controller  and  linear  mixer  in  each  channel.  Us¬ 
ing  the  gsdns  discovered  previously  through  pole  placement  for  the  PI  controller,  the  mixer 
gains  are  used  to  place  poles  that  yield  a  better  formation  response  to  lead  maneuvers. 


4-32 


■ 


These  linear  mixer  gains  are  used  in  the  nonlinear  formation  control  system  model  and 
are  refined  in  order  that  the  best  tracking  response  be  obtained.  These  newly  determined 
mixer  gains  are  then  kept  constant  while  the  PI  controller  gains  are  again  used  to  place 
the  poles  of  the  system.  The  gains  determined  through  this  experimentation  are  simu¬ 
lated  again  in  the  nonlinear  formation  control  model  and  refined  if  the  tracking  response 
is  inadequate  [2:page  5-12]. 

SISO  closed-loop  transfer  functions  are  obtained  for  the  y  and  x-channels  that  re¬ 
late  heading  and  velocity  to  their  respective  commands.  They  are  determined  using  the 
following  equation: 


G(s)  =  Ce,[s/-^,r'r„  (4.98) 

The  gain  adjustment  strategy  is  to  induce  the  fastest  poles  corresponding  with  the  largest 
residues  [2:page  5-13).  A  first-order  response  is  desired. 

Why  would  someone  want  a  first-order  response?  In  this  case,  the  true  system  is 
nonlinear  because  of  rate  limiters.  A  first-order  transfer  function  has  the  form  of 


Yi»)  = 

'  '  TS  -I- 1 

where  r  =  time  constant  (seconds). 

The  output  time  response  is 


(4.99) 


y(t)  =  1  -  e  r. 

Differentiating  the  output  time  response  yields, 

./  X  1-1 

y(t)  =  -«  '  • 


(4.100) 


(4.101) 


4-33 


The  system  in  this  thesis  is  nonlinear  due  to  the  inclusion  of  the  limiters  in  the  first- 
order  models.  However,  if  the  rate  values  in  the  model  never  reach  the  value  of  the  limiter, 
then  the  system  is  operating  in  the  linear  region  and  is  essentially  a  linear  system  as  long 
as  it  operates  in  this  region.  A  rate  is  equal  to  j .  Now,  if  a  rate  limiter  in  a  system  limits 
the  output  between  ±4.7°,  and  7  =  4,  the  system  wiU  be  linear.  Therefore,  it  is  in  the 
best  interest  of  the  designer  to  ensure  the  system  is  linear  by  establishing  a  dominant  pole 
that  allows  the  system  to  remain  in  the  linear  region. 


Figure  4.13  Comparison  of  Linear  and  Rate-Limited  Velocity  Responses  [2:Figure  5.2] 

Buzogany  discovered  in  his  thesis  that  the  “rate-limiters  in  the  aircraft  [channel] 
models  tend  to  produce  an  under-damped  oscillatory  response  for  high  integral  gain  values. 
This  response  is  caused  by  the  increase  in  tracking  error  induced  by  the  presence  of  the 
rate-limiters  [2:page  5-13].”  This  is  shown  in  Figure  4.13.  Since  integral  control  is  required 
for  tracking,  there  is  a  tradeoff  situation  existing  in  which  integral  control  is  needed  but 
too  much  wiU  cause  instability. 
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The  formation  parameters  and  gains  used  in  this  analysis  are  provided  in 
Ttble  4.3.  These  are  the  same  values  used  by  Busogany  in  his  analysis.  The 
PI  controller  gains  in  Table  4.3  are  the  gains  Busogany  determined  through 
pole  placement. 

IMPORTANT:  The  following  analysts  is  performed  using  dimensional  quantities 
(i.e.,  time  constants,  gains,  poles,  etc).  This  is  done  to  aid  in  the  comparison  of  the  results 
in  Busogany ’s  thesis  and  this  research. 

It  is  important  to  examine  how  the  gains  used  in  this  research  are  nondimensionalised. 
Using  Equations  (4.14)  and  (4.15),  the  following  gains  can  be  nondimensionalised: 


However,  the  gains  used  in  this  section  are  dimensional. 

4-3.  J  Linear  Mixer  Gain  Adjustment .  The  objective  of  this  analysis  is  to  deter¬ 
mine  a  set  of  linear  mixer  gains  for  both  the  y  and  x-channels  that  will  result  in  better 
formation  responses  than  are  possible  with  a  PI  controller  alone.  This  is  accomplished  by 
keeping  separation  error  signal  mixer  gains,  and  k^,  equal  to  one  (see  Sections  4. 1.1.1 
and  4.1.2. 1)  and  varying  the  remaining  mixer  gains,  k^  and  h,,  to  determine  if  a  dominant 
pole  can  be  found. 

4- 3.  J.  1  Y-Channel  Mixer  Gains.  The  Laplace  transform  of  the  investigated 
output  response  due  to  a  step  input  is 


—  *»p(0  (v.)  ~  (v.)  ~  (v.) 

=  K  (?)  *,  =  *,  =  1  *.  =  *.  (^)  4,  =  4,  =  1 


V’»'('»)  —  ~  — ~  I  , 

S  S  +  P2  S  +  P3 


S  -|-P4  S 


(4.103) 
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This  is  detennined  by  applying  Equation  (4.98).  Keeping  =  1  and  varying  results  in 
Figure  4.14.  Varying  does  not  result  in  a  dominant  pole.  However,  it  seems  k^  >  2  yields 
a  better  pole/residue  response.  Therefore,  this  is  the  starting  point  when  the  nonlinear 
formation  model  is  simulated. 


Figure  4.15  Linear  System:  X-Channel,  Vw  Velocity  Step  Response,  Poles  and  Residues 
Versus  k. 


4-3. 1.2  X-Channel  Mixer  Gaina.  The  Laplace  transform  of  the  investigated 

i 

output  response  due  to  a  step  input  is 


xr  /  X  1  .  ’■j  **3  *‘4 

Viv(s)  —  -  +  — ; -  H - - - 1 — - — 

a  a +P2  S+P3  <  +  P4 


(4.104) 
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This  is  determined  by  applying  Equation  (4.98).  Keeping  k,  =  1  and  varying  k,  results  in 
Figure  4.15.  It  seems  k,  >  I  yields  a  better  pole/residue  response.  Tberefcwe,  this  will  be 
the  starting  point  when  the  nonlinear  formation  model  is  simulated. 

4.3.2  Tried  and  Error:  Mixer  Gains  .  A  heuristic  process  using  the  starting  point 
provided  by  the  previous  section  is  initiated  with  the  nonlinear  simulation.  The  maneuvers 
used  to  exercise  the  system  are  a  10"  heading  change  and  a  25  ^  velocity  change  shown 
in  Figures  4.16  and  4.17,  respectively.  Keeping  k,  =  k^  =  I,  mixer  gains,  k^si  and  k,,  are 
varied  to  improve  formation  performance. 

Experimentation  with  k^  and  k,  yields  some  interesting  observations.  Varying  k^ 
has  no  apparent  effect  on  a  velocity  change,  k^  does,  however,  affect  a  heading  change. 
Observations  concerning  k^  reveal: 

•  Increasing  k^  allows  wing  velocity,  Vwt  uud  x  separation  to  peak  at  lower  values  than 
were  previously  obtainable  without  the  mixer.  However,  y  separation  is  adversely 
affected  by  increasing  k^  since  it  initially  decreases  before  matching  the  trajectory  of 
the  channel  without  a  mixer  in  place.  As  k^  is  increased,  the  y  separation  deviations 
decrease  with  greater  magnitude  before  coming  back  to  their  nominal  value. 

•  Wing  heading,  or  Hw  in  the  plot,  has  a  slightly  better  response  with  the  mixer 
in  place.  It  has  a  slightly  quicker  rise  time. 

Observations  concerning  k,  reveal: 

•  Increasing  k,  allows  z  separation  to  peak  at  lower  values  than  were  previously  ob¬ 
tainable  without  the  mixer  for  a  velocity  change. 

•  However,  as  k,  is  increased,  the  plot  of  z  separation  during  a  heading  change  peaks  at 
a  higher  value  than  before.  Therefore,  k,  has  to  be  chosen  to  yield  good  performance 
for  a  velocity  change  and  heading  change. 

The  preceding  observations  agree  mathematically  with  the  y  and  z-channel  decoupled 
state  equations.  Equations  (4.28)  and  (4.52).  The  y-channel  is  completely  decoupled  from 
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Figure  4.16  C-130A:  Diamond  Formation,  10  Degree  Heading  Change,  Time  Response 
Comparison  of  Buzogany’s  (solid  line)  and  New  Gains  (dashed  line) 
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the  c-chaonel.  Its  only  modelled  disturbance  is  the  commanded  lead  heading  change 
Therefore,  its  mixer  gains  (i.e.  in  this  case)  should  be  effective  only  for  a 
heading  disturbance.  However,  the  y-channel  is  coupled  into  the  s-channel  through  the 
disturbance  =  V’w,  -  ijfw  •  Hence,  heading  changes  cause  c  separation  and  wing  velocity 
disturbances  as  seen  in  Figure  4.16  and  shown  mathematically  in  Equation  (4.52),  the 
decoupled  z-channel  state  equation,  affects  the  coupling  into  the  c-channel,  or  more 
distinctly,  hy  affects  the  disturbance  into  the  c-channel.  Therefore,  must  not  be  too 
large  since  it  directly  impacts  the  magnitude  of  the  disturbance  in  the  c-channel. 


Figure  4.17  C-130A:  Diamond  Formation,  25  ^  Velocity  Change,  Time  Response  Com¬ 

parison  of  Busogany’s  (solid  line)  and  New  Guns  (dashed  line) 


Table  4.6  Diwimainnal  Mixer  Gains 


PARAMETER 

NEW 

PREVIOUS 

mmsmmi 

4 

2 

0 
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Table  4.7  DimensioualP 

Con, 

trQllet_Gauia_UBed  In  Simulatimu 

V 

.01 

1 

.0075 

.75 

In  the  same  manner,  k,  nuut  also  be  chosen  for  better  performance  in  response  to 
commands  in  both  channels.  It  has  to  ]rield  better  performance  during  a  velocity  change, 
while  allowing  the  c  separation  during  a  heading  change  to  be  small.  The  gains  used  for 
and  k,  are  listed  in  Table  4.6. 

It  can  now  be  concluded  that  the  addition  of  a  linear  mixer  improved  aircraft  per¬ 
formance. 

(NOTE:  It  is  an  interesting  observation  that  k^  only  has  an  effect  during  a  heading 
change,  while  k,  has  effects  both  during  heading  and  velocity  changes  due  to  the  coupling 
of  the  p-channel  into  the  x-channel.) 

4.3.S  Varying  Y-Channel  PI  Controtler  Gains  Only.  The  PI  controller  gains 
affecting  the  y-channel,  k^  and  k„,  are  now  varied  to  determine  if  a  dominant  pole  can 
be  found.  The  Laplace  transform  of  the  output  response  being  investigated  is 


7+ 


- 1 - ; — 

Pa  *  +  P3 


S  -I-  P4  s  4 


(4.105) 


As  in  Buzogany’s  thesis,  the  PI  controller  gains  are  varied  while  the  pole  and  residues 
are  determined.  The  mixer  gains  are  the  same  as  in  Ihble  4.6.  The  ratio  ^  is  fixed  while 
is  varied  from  10~^  to  10.  Figure  4.18  shows  how  varying  the  PI  controller  gains 
affects  the  poles  and  residues  of  the  output  response.  There  is  not  a  gain  range  that 
allows  a  dominant  pole.  The  problem  lies  with  including  the  lead  states  into  the  model. 
It  contributes  a  pole  at  -1.5  with  residue  at  approximately  1.  This  pole  interferes  with 
other  poles’  dominance.  Therefore,  varying  the  PI  controller  gains  alone  does  not  3rield  a 
dominant  pole. 
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Figure  4.18  Linear  System:  Y-Channel,  ^  Heading  Step  Response,  Poles  and  Residues 
Versus  Controller  Gains,  4“  =  100 
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Figure  4.19  Linear  System:  X-Channel,  Vw  Heading  Step  Response,  Poles  and  Residues 
Versus  Controller  Gains,  =  100 

••f 
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4-S.4  Vatying  X-Channel  PI  Controller  Gains  Only.  The  gains  affecting  the 
z-channel,  and  k^p,  are  now  varied  to  determine  if  a  dominant  pole  can  be  found.  The 
Laplace  transform  of  the  output  response  being  investigated  is 


lV(s)  =  i 

s  s+p. 


+ 


+ 


•+P4  »  + Pi 


(4.106) 


Like  the  y-channel,  the  PI  controller  gains  are  varied  while  the  poles  and  residues  of 
the  system  are  determined.  The  mixer  guns  are  those  used  in  Table  4.6.  The  ratio  ^  is 
fixed  while  k^p  is  varied  from  10~^  to  10.  Figure  4.19  shows  how  varying  the  PI  controller 
gains  affect  the  poles  and  zeros  of  the  output  response.  Once  again  there  is  not  a  gain 
range  that  yields  a  dominant  pole. 


Since  no  dominant  pole  can  be  found,  a  trial  and  error  process  is  used  to  determine 
the  best  formation  response.  Buzogany’s  PI  controller  gains  listed  in  l^ble  4.3  are  used  as 
a  baseline. 

4-3.5  Trial  and  Error:  PI  Contrxdler  Gains.  Alter  determining  the  mixer  gains 
in  the  previous  section,  attention  is  now  focused  on  the  PI  controller.  A  trial  and  er¬ 
ror  approach  is  used  once  again  to  investigate  the  possibility  of  improving  the  formation 
performance. 

Using  Buzogany’s  PI  controller  and  the  linear  mixer  determined  in  the  previous 
section  as  a  starting  point,  the  PI  controllo^  gains  in  both  channels  are  changed  while 
the  mixer  is  kept  constant.  The  gsdns  in  Table  4.8  yield  better  performance  than  that  of 
Section  4.3.2  when  implemented  in  the  simulation.  (See  Figures  4.20  and  4.21). 


Table  4.8  Adjusted  Dimensional  PI  Controller  Gains  Used  In  Simulations  To  Achieve 
Better  Performance 


kxp 

.1 

1 

.1 

1 
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Figure  4.20  C-130A;  Diamond  Formation,  10  Degree  Heading  Change,  Time  Response 
Comparison  of  Figure  4.16  (solid  line)  and  New  PI  Gains  (dashed  line) 
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Figure  4.21  C-130A:  Diamond  Formation,  25  ^  Velocity  Change,  Time  Response  Com¬ 

parison  of  Figure  4.17  (solid  line)  and  New  Gains  (dashed  line) 
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Similar  to  Section  4.3.2,  tht;  y-channel  gains  only  affect  a  heading  change  while  the 
z-channel  gains  affect  both  heading  and  velocity  changes.  Once  again,  this  is  due  to  the 
disturbance  in  the  e-channel  due  to  the  y-channel  coupling  into  the  c-channel. 

4.S.6  Comparison  of  New  DimensiontJ  Gains  With  Thoae  of  Dargan's.  The  gains 
determined  by  Dargan  were  found  empirically,  like  the  new  gains  listed  in  Tables  4.6  and 
4.7.  The  gains  Dargan  used  are  in  Table  4.9.  These  same  gains  were  used  by  Busogany 
in  Chapter  7  of  his  thesis  when  using  first-order  models.  C-130B  models  are  now  used  to 
compare  to  simulation  runs  accomplished  by  Busogany  in  Chapter  7  of  his  thesis. 


Thble  4.9  Dargan’s  Dimensional  Gains  (See  Also  Thble  4.10  For  Adjusted  k,) 


Dimensional  Gains 

Value 

.5 

.05 

.17 

.02 

*. 

2 

5 

1 

10 

As  can  be  seen  by  examination  of  Figures  4.22  and  4.23,  the  new  gains  provide  a 
better  response  for  all  states  of  the  system  with  the  notable  exception  of  Vw  when  a  30° 
heading  change  is  commanded  as  in  Figure  4.22.  The  gains  used  by  Dargan  and  Busogany 
have  a  much  better  response  for  Vw  than  the  new  gains  found.  The  question  to  ask  is  which 
provides  the  best  overall  performance.  The  constraint  most  important  to  formation  flying 
is  that  of  maintaining  an  adequate  separation  distance.  No  matter  what  the  maneuver  is,  if 
the  aircraft  involved  cannot  perform  it  in  a  timely  manner  due  to  their  formation  separation 
response  then  the  formation  flight  control  system  is  not  adequate.  The  new  gains  yield 
a  much  better  separation  response  for  both  heading  changes  and  velocity  changes.  The 
increase  in  wing  velocity  for  a  30°  heading  change  in  Figure  4.22  is  approximately  14 
^  more  than  the  response  of  the  Buzogany/Dargan  gains.  This  is  a  minimal  change. 
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Figure  4.22  C-130B:  Diamond  Formation,  30  Degree  Heading  Change,  Time  Response 
Comparison  Of  Dargan/Busogany  Gams  For  Mixer  And  PI  Controller  (solid 
line)  And  New  PI  Gains  (dashed  line).  First-Order  Models  Are  Used 
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Figure  4.23  C-130B:  Diamond  Formation,  25  ^  Velocity  Cliange,  Time  Response  Com¬ 
parison  Of  Dargan/Buzogany  Gains  For  Mixer  And  PI  Controller  (solid  line) 
And  New  PI  Gains  (dashed  line).  First-Order  Models  Are  Used 
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Examining  the  ]f-iq>aration  response  deviation  in  the  same  figure  shows  it  to  be  much 
better  (approximately  20  ft).  Sq>aration  between  aircraft  is  the  more  important  of  the  two 
indicators. 

Even  though  the  new  gains  give  a  better  overall  response,  it  seems  the  responses  can 
be  refined  slightly  more.  Therefore,  following  additional  trial  and  error  approach,  the  only 
gain  that  needs  to  be  changed  is  h,.  The  new  k,  is  listed  in  Table  4.10. 


Tnble  4.10  Adjttsted  Gain  To  Refine  Response  of  Figure  4.22 


Adjusted  Dimensional  Panmeter  Value 


vigure 

9 


Examining  Figure  4.24  it  is  seen  that  the  response  is  slightly  better  than  that  of 
Figure  4.22.  During  simulations  it  is  seen  that  there  is  a  tradeoff  between  Vw,  wing 
velocity,  and  c-separation.  As  Vw  increases  there  is  less  deviation  from  the  nominal  x- 
separation  distance  and  vice-versa.  This  makes  physical  sense.  If  a  diamond  formation 
composed  of  two  aircraft  makes  a  heading  change,  the  wing  aircraft  will  have  to  increase 
speed  if  it  is  on  the  outside  of  the  turn,  wing  speed  is  increased,  there  should  be  less 
deviation  in  the  nominal  z-separation  (following  this  same  argument,  the  converse  is  also 
true).  An  additional  benefit  is  the  operation  of  the  wing  velocity  stays  longer  in  the  linear 
region  of  operation  than  previously. 

4.4.  Dimensional  Eigenvalue  Comparison 

(See  Appendix  A  for  dimensional  plant  models)  Applying  Equation  (4.98)  to  the 
dimensionalized  y-channel  yields  the  y-channel  closed-loop  transfer  functions  which  when 
disturbed  by  the  lead  aircraft  heading  command  results  in  the  Laplace  output  response 
below. 


^(s)  =  G(s)V»x,e(<)  =  I  0  10  0  [s/-A,]  (4.107) 
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Figure  4.24  C-130B:  Diamond  Formation,  30  Degree  Heading  Change,  Time  Response 
Comparison  Of  Dargan/Busogany  Gains  For  Mixer  And  PI  Controller  (solid 
line)  And  New  PI  Gains  (dashed  line).  First-Order  Models  Are  Used,  k,  = 
3. 
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Figure  4.25  C-130B:  Diamond  Formation,  25  ^  Velocity  Change,  Time  Response  Com¬ 
parison  Of  Dargan/Buxogany  Gains  For  Mixer  And  PI  Controller  (solid  line) 
And  New  PI  Gains  (dashed  line).  First-Order  Models  Are  Used,  =  3. 
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where,  the  mixer  aad  PI  controller  gains  are  those  specified  in  Thble  4.11. 


Thble  4.11  Mixer  and  PI  Controller  Gains  Used  In  Dimensional  Eigenvalue  Comparisons 


Parameter 

k,i 

ksp 

kv 

Dimensional  Value 

.1 

1 

4 

1 

.1 

1 

3 

1 

Now, 


■  0 

0 

0 

375 

-1.5 

0 

375 

0 

-1.5 

-750  ■ 
1.5 

0 

.r,= 

0 

0 

1.5 

.1 

380.6 

369.4 

-756 

6 

(4.108) 


therefore. 


1  .00035587  -1.9689  .9685  .000013813 

■  s  s  +  .10003  s+  .74319  s  +  1.5  s  +  756.66 


(4.109) 


Busogany’s  Laplace  output  response  is  (see  Busogany  [2],  Equation  (5.58)) 


V^v(«)  =  G(s)t('i,(«)  = 


Busogany’s  Response 
1  .  .0001  .  -1.0010  .0010 
s  s  +  .Ol  s  +  .7492^ s  +  750.7408 

ilomtnant 


(4.110) 


As  can  be  seen  in  Equation  (4.109),  there  is  not  any  one  pole  that  is  dominant  as  in 
Equation  (4.110).  Therefore,  its  behavior  will  not  be  a  classical  first-order  response  when 
disturbed  by  a  step  input. 

An  interesting  observation  can  be  made  by  comparing  the  new  eigenvalues  and 
residues  (Equation(4.109))  to  those  of  Busogany  (Equation  (4.110)).  There  is  a  close 
correspondence  between  the  eigenvalues,  but  not  the  residues,  however.  The  eigenvalue 
that  does  not  correspond  to  Busogany’s  response  is  due  to  the  inclusion  of  the  lead  aircraft 
states  into  the  Formation  Flight  Control  problem  formulation.  Another  observation  is  that 
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since  the  sum  of  the  residues  equals  0,  the  denominator  is  two  or  more  degrees  higher  than 
the  numerator  [Srpage  117].  (This  last  statemeat  was  proven  to  be  true  using  the  ss2tf 
(state-space-2-transfer-functi<m)  command  in  Matlab) 

Applying  Equation  (4.98)  to  the  c-channel,  where 


A. 


■  0 

-1 

1 

0 

■  0 

750  ■ 

0 

-3 

0 

3 

r  — 

0 

0 

0 

0 

-3 

0 

3 

0 

.1 

7.7 

-7.7 

-9 

9 

750 

C.  =  [  0  1  0  0  ]  (4.111) 


results  in  the  following  Laplace  transform  of  the  output  response 


Vwis)  =  G(3) 
therefore,  with  Vx,,  as  the  input, 


VlA^) 

D^(s) 


=  [o  1  0  o]l*/-A.r'r.^  (4.112) 


,  1  .37130  -.63590  -1  .26469 

yw{»)  =  -  +  +  -  ;■  + 


s  s  +  .12299  8  +  .20906  s  +  3  s  +  11.668 

and,  with  as  the  input. 


(4.113) 


„  .  -  750  423.21  -1190.1  5.9348x-^<  16.862 

~  8  8  +  .12299  s  +  .20906  s  +  3  '^  8  +  11.668  (4114) 


Buzogany’s  Laplace  output  is  (see  Buzogany  [2],  Equation  (5.59)), 


Buzogany’s  Response 


T,  /  V  ^  /  X  rx  /  X  1  0138  -6.6677  5.6538 

Vw{»)  =  GA»)DA»)  =  -  -H  ■  ;,;,x,  +  -  -,-7:;  + 


8  S  +  .OIOI  s  + 1.3719  s  + 1.6179 


(4.115) 


4-54 


Once  again,  there  are  some  interesting  observations  that  can  be  made  by  comparing 
the  eigenvalues  and  residues  of  the  two  responses.  There  is  not  a  close  correspondence 
between  the  eigenvalues  of  the  different  approaches  as  there  is  in  the  y-channel.  There  is  no 
dominant  pole  in  the  new  response  and  Busogany’s  response.  The  residues  iq>prozimately 
add  up  to  sero,  thus  the  denominator  is  two  or  more  degrees  higher  than  the  numerator 
[5:page  117].  (This  last  statement  was  proven  to  be  true  using  the  ss2tf  (state-space-2- 
transfer-fimction)  command  in  Matlab).  Once  again,  the  first  order-model  pole  shows  up 
in  the  response.  This  pole  is  not  allowing  any  other  pole  to  become  dominant. 

The  disturbance  response  is  also  interesting.  It  has  a  large  magnitude  and  no  true 
dominant  pole.  Its  step  response  will  not  exactly  be  a  first  order  response,  but  will  be 
closer  than  the  velocity  and  heading  step  responses.  It  is  of  no  immediate  concern  to  the 
Formation  Flight  Control  system  since  it  is  the  response  of  the  disturbance,  V'We  ' 
which  goes  to  zero  in  steady-state. 

4-4- J  PloU  of  Linear  System  Response  To  A  Unit  Step.  The  linearised  model 
responses  are  simulated  in  order  to  compare  the  newly  derived  model  with  Buzogany’s 
model.  Figure  4.26  shows  the  ijfw  response  to  a  unit  step  disturbance  for  both  the  new 
system  and  Buzogany’s,  Figure  4.27  shows  the  Vw  response  to  a  unit  step  disturbance  for 
both  systems,  and  Figure  4.28  shows  the  Vw  response  to  the  y-channel  coupling  into  the 
x-channel. 

As  can  be  seen  from  these  figures,  the  addition  of  the  lead  aircraft  states  into  the 
linearized  plant  model  has  resulted  in  slower  rise  times  for  the  different  responses.  This 
makes  physical  sense.  With  addition  of  more  dynamics  into  the  system,  additional  phase 
lag  is  introduced,  thus  slowing  the  responses. 

4-5  Conclusion 

This  chapter  develops  a  linearized  and  parameterized  MIMO  plant.  It  then  presents 
a  thorough  investigation  of  x  and  y-channel  stability  and  their  stability  envelopes.  The 
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Figure  4.26  (Linear  System)  TOP:  New  ^  Response  to  a  Unit  Step  Disturbance  (^l.) 

with  Mixer  and  PI  Controller;  BOTTOM:  Busogany’s  ijfw  Response  to  a 
Unit  Step  Disturbance  with  Mixer  Only 
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Figure  4.27  (Linear  System)  TOP:  New  Vw  Response  to  a  Unit  Step  Disturbance 

with  Mixer  and  PI  Controller;  BOTTOM:  Buzogany’s  Vw  Response  to  a 
Unit  Step  Disturbance  (V^)  with  Mixer  Only 


Figure  4.28  Linear  System:  VW  Response  to  a  Unit  Step  Disturbance  {dip) 


4-58 


i 

I 


performance  of  the  new  system  and  previous  work  is  compared.  Finally,  an  eigenvalue 
comparison  between  present  and  previous  work  is  accomplished. 

The  final  and  overall  result  of  the  previous  analysis  is  that  the  addition  of  a  linear 
mixer  does  indeed  improve  system  performance. 
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V.  Three  Dimensional  Energy  Excursion  Minimizing  Maneuvers 


In  Chapter  IV,  formation  velocity  and  heading  changes  accomplished  in  the  plane 
(i.e.,  they  are  two  dimensional  maneuvers)  are  considered.  Consequently,  the  wing  air¬ 
craft  has  a  smaller  or  larger  turn  radius  than  the  lead  aircraft  during  a  heading  change, 
due  to  the  constraint  of  two  dimensional  maneuvers.  This  necessitates  a  wing  accelera¬ 
tion/deceleration  sequence.  This  chapter  investigates  the  possibility  of  performing  three 
dimensional  maneuvers  and  maintaining  formation  integrity  in  steady-state.  This  research 
is  a  continuation  of  the  research  begun  by  Buzogany  in  1992  [2:Chapter  6].  Busogany’s 
thesis  formulation  only  considered  automatic  flight  control  heading  change  maneuvers. 
This  research  is  extended  to  include  the  design  of  a  formation  hold  autopilot  using  energy 
nniniiniiiing  maneuvers  which  will  allow  velocity  and  altitude  change  maneuvers  in  addition 
to  heading  change  maneuvers. 

Three  Dimensional  Maneuvers  afford  the  possibility  of  conserving  the  wing  aircraft’s 
energy.  Instead  of  accelerating/decelerating  to  maintain  formation  integrity  as  in  the  two 
dimensional  problem,  the  wing  aircraft  now  has  the  ability  to  accelerate/decelerate  by  de¬ 
creasing/increasing  altitude,  respectively,  without,  however,  having  to  change  its  energy 
state.  The  latter  is  most  intimately  connected  to  throttle  activity,  which  is  reduced.  This 
is  done  at  the  expense  of  separation  requirements.  Therefore,  maintaining  the  vertical 
separation  requirement  between  the  lead  and  wing  aircraft  is  relaxed  in  order  to  accommo¬ 
date  maneuvers  outside  of  the  lead  aircraft’s  flight  plane.  Also  as  a  result,  the  h-channel 
is  coupled  into  the  x-channel. 

The  analysis  of  three  dimensional  maneuvers  is  accomplished  in  this  chapter  with 
first-order  aircraft  models. 
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5.1  Development 

Following  the  same  initial  theoretical  development  as  in  Biuogany,  Ch^>ter  6,  the 
(total)  specific  energy  of  an  object  is  the  sum  of  its  kinetic  and  potential  energy  per  unit 
mass.  Therefore,  the  total  specific  energy  for  the  wing  aircraft  is 

ew  =  (5-1) 

where 

g  =  acceleration  of  gravity 
ew  =  wing  aircraft’s  specific  energy 
=  wing  aircraft’s  altitude 

Differentiating  Equation  (5.1)  and  setting  iw  equal  to  zero  leads  to  the  fundamental 
relationship  for  energy  conservation, 

ew  =  VoVw+ghw  =  0  (5.2) 

Using  the  first-order  aircraft  altitude  response  model  from  Buzogany’s  thesis  [2:Eiquation 
(6.3)], 


hw  — - hw  H - Wcj  (5-3) 

and  substituting  Equation  (4.3),  the  first-ordo’  aircraft  velocity  response  model,  into  Equa¬ 
tion  (5.2)  yields  the  following  equation, 

ejy  =  — -  -^Viv  +  — Viv.,  (5.4) 

which  shows  the  total  specific  energy  development  when  it  is  controlled  by  an  altitude 
change  and/or  a  velocity  change. 
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Like  Chapter  IV,  it  ia  desirable  to  nondimensionalise  and  parameterise  Equation 
(5.4).  Once  again,  Equations  (4.14),  (4.15),  and  V„  will  be  used  to  nondimenaionalue 
Equations  (5.3)  and  (5.4)  (Equation  (4.3)  has  already  been  nondimenaionaliaed.  See  Equa¬ 
tion  (4.18)).  Therefore,  the  additional  nondimenaional  variables  are: 


W  =  ‘f'ktr 

«iv  =  Civ  =  ^ 


(5.5) 


(The  symbol  u  now  dropped  for  convenience.  All  variables  are  nondimen^ 
sional  unless  otherwise  noted.) 

The  nondimenaional  form  of  Equation  5.4  is  , 


A  common  factor  in  Equation  (5.6)  is  This  expression  is  nondimenaional.  Therefore, 

(51) 

where, 

(  =  nondimensional  altitude  to  velocity  proportionality  constant  [2:page  6-2]. 
Substituting  (  for  M  yields. 


ew  =  - hwr  -I - Vw  +  ^ -I - (5.8) 

“^hw  ^V,r  ’*'aw 

With  a  slight  abuse  of  notation,  the  derivative  states  of  Equation  (5.2),  can  be 
assumed  to  be  perturbations,  and  Vo  is  the  unperturbed  formation  speed.  Equation  (5.2) 
now  becomes 
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ew  =  gfiw  +  VoVw 
dimensional  equation 

Rearranging  Equation  (5.9)  and  isolating  Vw  yields, 


Vw  = 

dimensional  equation 

In  nondimenaional  form, 


(5.10) 


Vw  = 


(5.11) 


Now,  Equation  (5.11)  is  substituted  into  Equation  (5.8).  The  resulting  equation  is 


ew  = - ew  +  ^  ( - J  hwr  +  ^ - hw,  + - Vw, 

‘’'Vw  \‘fVw  ^ir/  ‘^w 


(5.12) 


Taking  Equation  (4.17),  the  nondimensionalised  i  separatimi  equation  of  Chapter  IV,  and 
substituting  Equation  (5.11)  for  Vw  results  in  the  following: 


*  =  (hw  ~ 


(5.13) 


where. 


,  „  sma,,  ,  , 

d  =  Vji  +  -  (^^e  ~  ^ ) 


In  addition,  the  lead  aircraft  altitude  state  is  avsulable: 


•  11. 
hi  =  4-  — 


(5.14) 


(5.15) 
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Therrfove,  h,,  the  altitude  error  Bigual  can  be  generated. 

Equations  (5.3),  (5.12),  (5.13),  and  (5.15)  can  now  be  put  into  state  space  form, 


4"  r„^Xl 

hence. 


X 

o 

1 

V/ 

o 

1 _ 

X 

hw 

o 

1 

o 

o 

hw 

ew 

0  - *-)  — L  0 

'  \^w 

ew 

.  . 

0  0  0 

”  L. 

■  0  o' 

1  0 

J 

0  ^ 

■  Vw, 
hw 

+ 

0  0  0 

0  0  0 

■  Vl  ■ 
hL, 

0  0 

0^0 

L  J 

. 

where,  as  in  Chapter  IV,  Section  4.1.2. 


(6.W) 


(5.17) 


As  in  Busogany’s  thesis,  there  is  a  need  for  a  PI  controller  to  be  designed  to  com¬ 
mand  the  wing’s  altitude  response  to  a  change  in  x-separation  [2:page  6*2].  Hence,  the  PI 
controller  is 


(5.18) 


where, 

=  nondimensional,  three  dimensional  x-separation  error,  proportional  gain 
^sst  =  nondimensional,  three  dimensional  x-separation  error,  integral  gain 
Differentiating  Equation  (5.18)  results  in 


—  hxj^X  hsxiX 


(6.19) 
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Sttbatitttting  Equation  (5.13)  into  Equation  (5.18)  yiddf  a  result  limilar  to  Butogany 
[2:Equation  (6.10)], 


hw,  =  +  kttiZ  (5.20) 

Now  a  new  control  concept  it  needed  to  relate  energy  contervation  and  commanded 
velocities.  This  new  control  concept  it  baaed  on  using  Vw,  as  a  control  to  conserve  the 
wing’s  total  specific  energy  during  a  formation  maneuver.  Conserving  the  wing  aircraft’s 
total  specific  energy  implies  moderate  throttle  usage.  Therefore, 


ew 


V{T  -  D) 


where, 

T=  aircraft’s  thrust 
D  =  urcraft’t  drag 
m  =  aircraft’s  mass 
Using  a  PI  control  law  as  before. 


(5.21) 


Vw^  =  kstpvf^  +  hx*»v  /  he  (ft  (5.22) 

Jo 

where, 

h,  =  hi  -  hw  =  altitude  error 

kg,py  =  nondimensional,  three  dimensional  altitude  error,  proportional  gain 

htMiv  =  nondimensional,  three  dimensional  altitude  error,  integral  gain 

However,  Equation  (5.22)  does  not  conserve  the  wing  aircraft’s  energy.  If  iw  i>  *«t  equal 
to  sero  in  Equation  (5.12),  the  equation  is  solved  for  Viv.,  and  the  new  relation  for  Vw,  is 
added  to  Equation  (5.22),  the  following  control  law  is  obtained: 
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VW.  =  fVw  { - i - fcw  -  i +  *»*«v  f 

\TVw  I'^iir  A'wJ  ^ir  /  .  Jo  _ 

>  I  I  ■■■  '■  V  -  ^  ^  ■"'  '  V  "— ^ 

Energy  Conserviiig  =  elevator  usage  Energy  Changing  =  throttle  usage 

(5.23) 

This  control  law  is  a  contradiction  in  terms,  because  there  are  two  distinct  and  different 
parts.  The  part  determined  by  setting  ew  =  0  conserves  the  wing  aircraft  specific  energy, 
and  thus  works  the  elevator  of  the  wing  aircraft.  In  contrast,  the  PI  controller,  contingent 
upon  altitude  error,  changes  the  wing  specific  energy.  It  is  actually  commanding  the 
throttle  usage  of  the  wing  aircraft.  The  two  pieces  form  an  energy  minimiihig  control  law. 
Differentiating  Equation  (5.23)  yields: 


Vw,  =  ( - iw  +  i  [ - 1  hw  -  f - Aht.)  +  (5-24) 

\Tv,r  LTlktr  "^tr  J  ^kw  / 

The  closed  loop  xs-channel  is  found  by  augmenting  the  open  loop  state  space.  Equa¬ 
tion  (5.17),  with  Vw,  and  hw,-  The  resiilting  closed  loop  state  space  is: 


~  "I"  (5.25) 

n.d  =  (5.26) 


hence. 
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X 

hw 

ew 

ht 

Vw, 

where, 


0 

-1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

’^ir 

0 

'  V*Vw  J 

_ i_ 

0 

1 

0 

0 

0 

1 

0 

0 

Du. 

1 

Dhw, 

h*»< 
r  1 

h,,pi 

n 

0 

0 

1 

0 

+ 


0 

0 

0 


0 

0 

1 
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0 

0 


•mnp 


0 


nma 

^,r. 

L  tima 


n.  = 


1  0  0  0  0  0 
0-10100 


X 

K 


Vt 

^Lc 


X 

hw 

«w 

hi, 

Vw, 

.  hw.  . 

(5.27) 

(5.28) 


V^Vir  "^hw  /  \^k»r 


~^^9Mpv  hjtjjy. 

*^hw 

1  .  1  . 

+  ^«rfz — *«»r 

^wr  ^hw 


—  —  -4-  h 


>1  ^  W  1  Ml.  1 

~  1-  ~  IV  ) 

^fc»r  ^*w  '^fctr  ^w '  ^hw 


(5.29) 

(5.30) 

(5.31) 

(5.32) 


and, 


=  augmented  zz-channel  state  vector 
=  augmented  xz-channel  plant  matrix 
Fgt^  =  augmented  «z-channel  input  matrix 
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=  augmented  cz-diamiel  input  vect<v 
^■>«i  =  augmented  sz-channel  output  vector 
Cat^  =  augmented  x2-channel  output  matrix 

5.1.1  Difference  Between  The  New  Approach  and  Prenoae  (Bnzogany’a)  Approach. 

Loi  Chapter  6  of  Buaogany’a  theaia,  eneigy  conierving  maneuvers  were  considered.  They 
were  accomplished  by  uning  the  control  hw.  end  feeding  back  z  separation.  In  the  new 
formulation,  energy  excursion  minimiaaticaa  is  investigated,  hw,  u  still  a  control  and  x- 
separation  is  fed  back,  however,  the  altitude  error,  is  also  fed  back  and  an  additional 
control,  Vtir.,  is  implemented. 

Previously,  during  a  heading  change,  wing  aircraft  velocity  changes  were  made  by 
changing  altitude  during  a  turn.  Once  the  turn  was  completed  the  wing  aircraft  used  its 
kinetic  energy  to  return  to  the  lead’s  flight  plane.  However,  if  the  lead  aircraft  actually 
changed  speed,  then  the  wing  aircraft  would  increase  or  decrease  its  own  velocity  by  chang¬ 
ing  its  altitude  and  thus  maintaining  constant  specific  enogy,  ew  Therefcne,  the  previous 
system  was  satisfactory  during  formation  heading  changes  but  not  for  formation  energy 
changes,  where  changing  velocity  or  altitude  meant  a  corresponding  change  in  the  other. 

The  concept  of  energy  excursion  minimising,  presented  in  this  duqyter,  tries  to  con¬ 
serve  energy  but  also  allows  the  wing  aircraft  to  track  the  lead  aircraft  through  both 
heading  and  velocity  changes.  This  is  accomplished  by  keeping  Buxogany’s  approach  in 
place  and  augmenting  it  with  a  PI  control  law.  Equation  (5.23)  without  the  PI  controller 
is  iw  =  0.  Therefore,  by  including  the  PI  controller  in  Eiquation  (5.23)  the  requirement 
for  keeping  the  iw  =  0  has  been  relaxed.  No  longer  does  energy  conservation  apply.  The 
object  of  the  xz-control  system  is  to  minimise  the  energy  swings,  namely  throttle  usage. 
Therefore,  during  a  heading  change,  the  aircraft  will  still  alter  its  x-separation  by  changing 
its  altitude,  but  it  will  return  back  to  its  flight  plane  quicker  due  to  moderate  throttle  us¬ 
age.  However,  energy  changes  can  now  be  accommodated.  Hence,  the  energy  minimising 
is  more  capable  than  energy  conserving  maneuvers.  However,  it  is  not  as  energy  efficient 
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during  formation  maneuvers.  Thus,  wing  aircraft  capability  has  been  increased  at  the 
expense  of  wing  aircraft  energy. 

5. 1.1.1  Why  The  Addition  of  The  PI  Controller  Works.  Looking  at  Equa¬ 
tions  (5.18)  and  (5.23),  one  might  wonder  how  a  control  scheme  can  be  effective  in  which 
one  control  law  commands  an  altitude  change  in  response  to  x-separation,  while  another 
control  law  commands  a  velocity  change  baaed  upon  altitude  error.  This  is  accomplished 
due  to  the  time  constants  of  the  two  controlled  states.  The  velocity  time  constant  is  smaller 
than  the  altitude  time  constant.  Therefore,  the  velocity  response  will  be  quicker  respond¬ 
ing  to  a  command  than  the  altitude  response.  Examining  Equations  (3.1)  and  (3.3),  and 
inserting  their  respective  time  constants  yields  the  following  transfer  functions: 

First-order  velocity  transfer  function  and  time  response  equation, 

V(.)  =  (5.3J) 

t>(f)  =  3«-’*  (5.34) 

First-order  altitude  transfer  function  and  time  response  equation. 


BW  =  (5.36) 

h(t)  =  .5e-“  (5.36) 

By  looking  at  Equations  (5.33)  and  (5.35)  or  Equations  (5.34)  and  (5.36),  it  is  obvious 
that  the  velocity  equations  will  reach  steady-state  faster  than  the  altitude  equations.  To 
illustrate  this  further,  look  at  Figure  5.1. 

Theoretically,  the  sequence  of  events  for  a  velocity  change  should  appear  in  this  order: 
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1 


Figure  5.1  Comparison  of  v(t)  (solid  line)  and  h(t)  (dotted  line) 

1.  During  a  maneuver,  the  x  separation  response  will  deviate  from  its  nominal  value 

2.  An  altitude  command  based  on  the  z-separation  error  will  be  issued.  This  command 
also  is  included  with  wing  altitude  perturbation,  hw,  to  command  the  velocity  chan¬ 
nel  first  order  model  to  increase  speed  m  order  that  the  z-separation  error  will  return 
to  sero.  This  is  the  energy  conserving  part  of  the  control  law 

3.  In  addition  to  this,  a  PI  controller  is  included  in  the  wing  velocity  command  control 
law  that  is  active  once  there  is  an  altitude  error  between  the  lead  and  wing  aircraft. 
Therefore,  in  addition  to  the  energy  conserving  part  of  the  control  law,  there  is  a 
energy  changing  part  which  commands  more  velocity  as  the  altitude  error  increases 

4.  Because  of  the  time  constants  of  the  two  channels,  the  altitude  response  should 
generally  be  slower  than  the  velocity  response  to  the  respective  commands.  Therefore, 
one  can  expect  an  initial  change  in  altitude  but  less  altitude  loss  due  to  the  increase 
of  velocity  during  an  energy  minimising  maneuver 

This  is  validated  by  simulation. 
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5.1.2  Static  Stability  Analysis .  As  in  Chapter  IV,  Section  4.1.2,  the  linear  system 
is  investigated  to  determine  if  the  system  is  able  to  track  a  step  disturbance  input  in  Vi  or 
hi^  (Remember  — »  0  at  steady  state.  See  Section  4.1.2).  It  is  imperative  that 

the  aircraft  be  able  to  track  an  altitude  or  velocity  command.  Like  the  two  dimensional 
case,  this  requirement  csm  be  stated  as  follows: 

for  disturbance  hi,^{s)  =  j. 


As  t  — ►  oo,  *(oo)  =  Vw^  =  0,  and  hw{oo)  =  ^^(oo)  =  hw^  -  f?rj  =  1, 

and  e(v(oo)  =  { 


(5.37) 


for  disturbance  =  7, 


As  t  — ►  00,  *(00)  =  hiv(oo)  =  hi,(oo)  =  hiv.(oo)  =  0, 
and  eHr(oo)  =  VUr.(oo)  =  ^1,(00)  =  1; 


(5.38) 


The  requirements  of  Equations  (5.37)  and  (5.38)  are  verified  by  observing  the  asz- 
channel  response  to  unit  step  distvirbances,  and  d^  (the  mathematics  in  this  section 

are  accomplished  with  Mathematica  [7]): 


x{t)„  =  lims-(5/-  A„)  ^r„  = 


■  0  0 
0  1 
1  ^ 
0  1 
1  0 
0  1 


0 

0 

MBg 

0 

■ino 

0 


(5.39) 


Since  — >  0  in  steady  state,  the  third  column  of  Equation  (5.39)  is  of  no  consequence 
and  is  disregarded  in  emalysis.  However,  the  first  two  columns  are  the  expected  results  from 
the  conditions  of  Equations  (5.37)  and  (5.38).  The  results  of  Equation  (5.39)  are  true  if 
Ag,  is  invertible.  This  implies  the  following  relation: 
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*.«  /  0  (5.40) 

Conclusion:  Once  again  (see  Chapter  IV,  Sections  4.1.2  and  4.1.1),  integral  control 
is  necessary  for  this  technique  to  work.  Furthmnore,  x  feedback  is  necessary. 

5.1.3  Dynamic  Siabiiity  Analysis  .  The  characteristic  equation  is 

det(5/  -  )  =  s®  +  +  js*  +  ks®  +  la*  +  ms  +  n  (5.41) 

where  the  coefficients  are. 
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-L  + 
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^Sw 

(5^  ■ 
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(  kggpkjixpv 

\  Tv^ 

1 

(  kgfiy  kxMf 
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1  1 

{  kxMfkxMiv 

kxMii 

kxxiki 


Tfc^  / 


•TMt"'XMpv 

k  ■  k 
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^»»j»kx»pv 

kxMjkxMpv  ^ 
’’’hL’TVyr  ) 


T-fc.  / 


kxMikx 


(5.42) 

(5.43) 

(5.44) 

(5.45) 

(5.46) 

(5.47) 


The  Routhian  array  for  Equation  (5.41)  is  in  Table  5.1. 

The  following  relations  result  from  the  Routh  criterion  application  of  Table  5.1: 
row  a® 


Tfcj  +  Tik„.  >  0 


(5.48) 
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Table  5.1  XZ-Channel  BxtuthiaD.  Array 


•i* -**-»*!+<«• 


+  ■nki.TVH-  +  >  0  (5.49) 


~{^SMp^wMpv'^ki,‘^kw)  ^*Miv^**p‘^ki^w  '1'  ^*n^»»fyr'^ki,'^w  ^PMp^Ktpv'^ht'^w'^w 
+fc»*tv^x»p^*t'^fc»r'A^»r  "I"  I^XMi^SMpv'^ht'^kw'^w  ~  +  ^**p^s*pv‘^hir'^w 

—  4x»»4«»tV^fct'^fct|r^W  "f"  ^XMp^**rV^hl,'^w(  ~  ^XMp^**Pv^ki‘^kw'^w( 
"■24*,jAxxptx»py.T(^^Tfc^7V^^  ~  ^xtp^kt/^Vw^  ~  ^**i‘^kt,‘^kw'^Vw^ 

+k,MiTjl„T^„^  +  -  k,,ikx,pTljhwrl^i^ 

~k\tiT^^Tu^Ty^^  >  0  (5.50) 

rows  are  in  Appenduc  B 


kxxikxxiv  ^  ® 


(5.51) 
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Equations  (5.48)  -  (5.51)  were  found  in  part  through  the  use  of  Mathematica  [7].  The 
equations  for  rows  are  very  large  and  difficult  to  analyse.  Due  to  time  constraints, 

the  stability  envelopes  of  these  stability  relations  will  net  be  explored  in  this  thesis  effort. 

These  stability  relations  show  the  interrelationships  between  the  controller  and  mixer 
gains,  first-order  aircraft  model  time  constants,  and  stability.  It  is  interesting  to  note  that 
this  analysis  agrees  with  the  static  stability  analysis  and  demonstrates  integral  control 
action  is  once  again  required  for  stability.  Equation  5.51  shows  that  kg,i  /  0  and  k^tiy  0. 
Therefore,  it  is  a  requirement  for  Energy  Excursion  Minimizing  Maneuvers  to  have  integral 
control  action  in  both  the  and  the  Vw^  control  laws. 

5.2  MathematiaU  Analyaia  of  Three  Ditneruional  Maneuvers 

Buzogany’s  thesis  effort  did  not  calculate  Vyf^ .  Instead,  V-w  was  determined  directly 
through  the  physics  of  the  problem  (i.e,  Vw  was  found  by  assuming  =  0  and  solving  for 
Vw  in  Equation  (5.11)).  It  is  useful  to  calculate  Buzogany’s  actual  velocity  control  law. 

Beginning  with  the  first-order  equation  for  the  aircraft  velocity  [Equation  (4.18] 

Kw  =  - — Viv  +  —  Vht.  (5.52) 

Substitute  the  relation  for  IV,  when  e^r  =  0  from  Equation  (5.11), 

Vw  =  —$hw  (5.53) 

Equation  (5.52)  becomes, 

-  i'hw  =  i—hw  +  — VW.  (5.54) 

Substitute  the  nondimensional  form  of  Equation  (5.3)  into  Equation  (5.54)  and  solve  for 
IV.  to  obtain 
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(5.55) 


Vw,  =  ihiw  ~  ~ 

\rkw  /  ^ir 

This  differs  from  Equation  (5.23)  in  two  ways: 

1.  The  PI  controller  utilising  altitude  error  is  absent.  This  was  not  in  Busogany’s 
formulation  of  the  problem 

2.  (ew  term  is  also  absent.  This  is  because  Busogany  assumed  ew  =  0 
Thus, 

Preposition  V.l  In  the  linear  caae,  the  new  control  law  in  Equation  (S.2S)  with  the  PI 
controller  gains  set  to  zero  will  be  the  same  velocity  control  law  Buzogany  had  previoudy 
in  his  thesis. 

However,  there  is  a  problem  with  Equation  (5.55).  This  was  Buzogany’s  Vw,  if  he 
commanded  the  first-order  velocity  model  given  in  Equation  (5,52),  However,  he  did  not 
do  that.  Instead,  he  used  Equation  (5.5S)  to  find  velocity.  In  order  to  determine  wftaf 
his  actual  Vw,  was.  Equation  (5.55)  needs  to  be  reexamined  and  all  terms  containing 
should  be  ne^ected.  This  leaves  the  following  equation: 

Vw.,  =  —^hw  (5.56) 

Which  is  the  same  as  Equation  (5.53).  This  equation  u  Busogsuiy's  true  Vw,  since 
he  did  not  have  the  first-order  velocity  model  included  in  his  simulation. 

The  differences  in  implementation  between  Busogany’s  and  the  new  simulation  are 
illustrated  in  Figures  5.2  and  5.3,  respectively. 
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Figure  5.2  Energy  Conservation:  Busogany’s  Three  F  'snsional  Maneuvers  Simulation 
Formulation,  With  the  Aircraft  Specific  Energy  Rate  =  0 


*•= 


Figure  5.3  Energy  Excursion  Minimisation;  New  Three  Dimensional  Maneuvers  Simula¬ 
tion  Formulation 
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5.S  GmphicttI  Comparison  of  Implementations:  Energp  Excursion  Minimixing  and  En¬ 
ergy  Conserving 

According  to  Section  5.2,  the  only  difference  between  energy  excursion  fniwimiaifig 
and  conserving  is  the  addition  of  the  PI  controller  and  the  tw  term  into  the  Vw,-  This  is 
verifiable  by  "turning  the  PI  controller  off”,  i.e.  set  PI  gains  equal  to  sero,  and  simulating 
energy  conserving  and  energy  Tninitnising  maneuver  control  techniques.  The  implementa¬ 
tions  are  depicted  in  Figures  5.2  and  5.3.  The  simulations  are  shown  in  Figure  5.4.  As  can 
be  seen,  the  time  responses  do  not  match.  Evidently,  the  nonlinear  simulation’s  limiters 
are  interfering  with  the  responses.  As  a  check,  the  interfering  limiters  are  eliminated  and 
the  resulting  responses  are  inspected.  Looking  at  Figure  5.5,  the  simulations  do  agree 
when  they  remain  within  the  linear  region  of  operation.  Therefore,  Preposition  V.l  has 
been  proven  correct  experimentally. 

This  comparison  of  implementations  shows  that  energy  conservation  is  really  just 
a  special  case  of  energy  minimising  with  the  PI  controller  gains  set  equal  to  sero.  The 
remainder  of  this  chapter  compares  the  two  methods.  In  the  remainder  of  this  thesis,  the 
term  “energy  conservation”  denotes  energy  minimisation  with  the  PI  controller  gains  set 
equal  to  zero. 

5.4  Energy  Excursion  Minimizing  Maneuvers  Simulation 

Once  agiun  the  simulations  are  accomplished  with  dimensional  gains.  Therefore,  it 
is  useful  to  determine  how  the  nondimensional  PI  controller  gains  in  Equation  (5.23)  are 
related  to  their  dimensional  counterparts. 

The  notation  is  dropped  in  the  rest  of  this  chapter.  The  reader  will 
be  told  when  equations  are  dimensioned  or  nondimensionsd.  This  section  deeds 
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ikw(ft)  dl(^ 


9 


Figure  5.5  Nonlinear  Simulation^  C-ISOA  Aircraft,  With  Velocity  and  Altitude 
Error  Signal  Limiters  Eliminated:  Comparison  of  Busogany’s  Implemen¬ 
tation  (solid  line)  and  New  Implementation  (dashed  line)  in  Equation  (5.23) 
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with  dimenaional  simiilalion  ruiia.  Thm^bre,  the  geiiis  are  ctwieiiwoeef  in  thb 
section. 

The  initial  conditioiu  for  the  three  dimen tiond,  energy  ctmserring  manenvera  ia  the 
lame  ai  Ihble  4.3.  The  PI  controller  gaini  for  Vw,  ***  in  their  respective  sectunu. 

5.5  Dimensional  Eigenvaine  Analy$u 

(See  Appendix  A  for  dtmeniuma/ plant  modeli)  Uiing  Bnsogany’i  results  as  a  starting 
point,  only  the  PI  controller  gains  associated  with  Vir.,  and  k,tiy,  are  varied.  The 
objective  of  this  section  is  to  compare  the  new  Lq;>lace  Transform  of  the  step  response  with 
Busogany’s  response  [2:page  6-5].  In  order  to  do  this,  the  output  matrix,  C,,  has  been 
been  changed  to  reflect  Busogany’s  output,  and  thus  examine  the  wing  vdocity  response. 
Therefore,  the  nondimenaional  wing  velocity’s  Laplace  output  step  response  is: 

Vl,r(s)  =  [0  -  o  0  0  0)  [•!  -  A..]-*  r„  -  (5.58) 

>  - - ^ 

Nondimensicmal 

Correspondingly,  the  dimensional  wing  velocity’s  step  response  is: 


Vw{») 


Vw{») 


0  -  Vb  0  0  0 
^0 


[s/-A„r^r..i 


Dimensional 
r2  ,  Ti  .  U 


»  «  +  Pa  *+P3 


+ 


+ 


re 


S  +  P4  »  + Pi  »  +  Pt  »  +  Pr 


(5.59) 

(5.60) 


As  in  the  previous  chapter,  it  is  desired  to  perform  an  eigenvalue  analysis  to  determine 
if  a  fast  response  is  possible  for  Vw  This  is  accomplished  experimentally  by  keeping  the 
ratio  of  =  100  while  varying  their  respective  values.  Keeping  this  ratio  equal  to 
100  is  an  engineering  “rule  of  thumb.”  Keeping  the  ratio  equal  to  100  is  not  a  proven 
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mathematical  concept.  However,  this  has  demonstrated  good  performance  in  the  past. 
Therefore,  it  is  used  in  this  analysis. 

Examining  Figure  5.6  shows  that  picking  gains  at  k„py  =  .0005  and  kg,iy  =  .000005 
yields  a  response  where  the  poles  p4  and  p5  are  dominant  (they  are  complex  coiyugates 
of  each  other).  These  gains  allow  for  tracking  of  a  step  input,  dominance  of  a  complex 
coiyugate  pair,  and  stability,  <  .005  for  stability  (See  poles  p4  and  p5). 

Previously,  Buxogany  found  the  wing  velocity  response  to  be: 


Dominant  Poles 

1  .043  '.72592Z  - 135.92‘  .72592Z135.92* ' 

s  s  +  0.0104  s  +  .2448  -  >.252  '''  s  +  .2448  +  >.252 

Busogany’s  Response 


(5.61) 


(7%ts  is  different  from  Equation  (6.Z1)  page  6-5  of  Bmogany’e  thesis.  The  etpreasion  in 
that  thesis  is  incorrect.  The  residues  do  not  add  up  to  zero,  which  is  a  requirement  for 
a  partial  fraction  expansion  whose  denominator  is  S  degrees  higher  than  the  numerator 
[5:page  117J). 

Transfer  function  and  Partial  FVaction  expansion  of  the  step  response  with  kg,fy  =  .0005 
and  ksxtv  =  .000005  yields: 

.1287(s  +  .5)(s  +  .00874  ±  >.00992)(s  +  .01) 

(s  +  .5)(s  +  .235  ±  >.243)(s  +  .00942  ±  >.00997)(s  +  .0104)  ^  ^ 

1  .050670  .040758 ->.040545  .040758 -|- >.040545 

s  s  +  .0104  s  +  .00942  +  >.00997  s  +  .00942  -  >.00997 

-.56609 ->.54353  -.56609  +  >  54353  0  _ 

s  +  . 235 +  >.243  s+. 235 ->.243  s  +  .5  ^  ’ 

''  V - ' 

Dominant  Poles 

These  gains  can  now  be  utilised  by  the  nonlinear  simulation  to  determine  if  they 
produce  desirable  results. 


VW(s) 

Vi(5) 

Viv(s)  = 
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5.5.1  Nonlinear  Three  Dimensional  Maiunvers  Simulation.  The  gains  used  for 
the  full  nonlinear  simulation  are  in  Table  5.2.  These  were  determined  by  the  eigenvalue 
analysis  in  the  previous  section.  Included  are  the  values  used  by  the  other  gains  in  the 
simulation. 


Table  5.2  PI  Controller  Values  For  Equation  (5.23)  and  Nonlinear  Simulation  Shown  in 
Figures  5.7  -  5.9 


Parameter 

Dimensional  Value 

^*MpV 

.0005 

k  . 

*SJ«V 

.000005 

1 

kpi 

.01 

hpp 

.75 

.0075 

1 

0 

k. 

1 

kp 

0 

The  vidues  from  Table  5.2  are  implemented  in  the  nonlinear  simulation  resulting 
in  Figures  5.7  -  5.9.  Examining  Figure  5.9,  where  a  25  ^  velocity  command  has  been 
issued,  shows  the  reason  for  adding  the  PI  controller  to  the  energy  conservation  part  of 
the  velocity  control  law.  Keeping  ew  constant  causes  the  wing  aircraft  in  the  energy 
conservation  response  to  decrease  altitude  wh«i  a  velocity  change  has  been  commanded. 
Allowing  ew  to  fluctuate  causes  the  wing  aircraft  utilising  energy  minimisation  to  recover 
back  to  its  nominal  altitude.  The  PI  controller  is  commanding  the  throttle  and  changing 
the  wing  aircraft  energy  state.  This  results  in  less  deviation  in  x-separation  than  the 
energy  conservation  approach.  In  Figure  5.9,  the  specific  energy  of  the  wing  aircraft  for 
both  methods  has  been  plotted  as  a  function  of  time. 

The  heading  changes  in  Figures  5.7  and  5.8  show  that  some  performance  is  traded 
for  more  capability.  The  rate  limiters  impact  the  response  here  dramatically.  The  energy 
conservation  responses  are  better  than  the  energy  minimizing  responses. 
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Figure  5.7  C-130A:  Diamond  Formation,  lO"  Side  Step  Heading  Change.  Nonlinear 
Simulation:  Energy  Conserving  Implementation  (solid  line)  and  Energy  Ex¬ 
cursion  Minimising  Implementation  (dashed  line) 


iinv($  iU(^  ,  H«(deg)  HKdeg) 


time 


Figure  5.9  C-130A:  Diamond  Formation,  25  ^  Velocity  Change.  Nonlinear  Simula¬ 
tion:  Energy  Conserving  Implementation  (solid  line)  and  Energy  Excursion 
Minimizing  Implementation  (dashed  line) 
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As  seen  in  Figures  5.7  -  5.9,  allowing  the  specific  energy  of  the  wing  aircraft  to  make 
excursions  away  from  its  nominal  value  through  the  addition  of  a  PI  controller  in  the  Vw. 
command  yields  two  conflicting  results  when  no  other  parameters  are  changed: 

•  Heading  changes  have  worse  responses  (i.e.,  there  is  undershoot  in  the  x  separation 
response  that  results  in  decreasing  aircraft  separation  ) 

•  Velocity  changes  have  much  better  responses  (i.e.,  x  separation  response  away  from 
nominal  is  less  and  the  aircraft  are  at  the  same  altitude  in  steady-  state  unlike  energy 
conserving  maneuvers 

5.6  New  Approach  To  Viv, 

Theoretically,  energy  excursion  minimization  should  result  in  more  aggressive  wing 
aircraft  response  to  lead  maneuvers  and  better  performance  than  energy  conserving  ma¬ 
neuvers.  A  better  response  can  be  found  in  velocity  changes  for  the  new  formulation, 
however,  there  is  a  worse  response  for  heading  changes.  Obviously,  the  present  formula¬ 
tion  for  energy  minimizing  maneuvers  is  not  as  good  as  was  previously  hoped.  The  concept 
needs  to  be  reexamined. 

The  whole  premise  of  energy  conserving  maneuvers  is  to  have  zero  rate  of  change  in 
the  perturbation  of  the  specific  energy  of  the  wing  aircraft.  Energy  excursion  minimizing 
maneuvers  allow  the  wing  aircraft  specific  energy  perturbation  to  make  excursions  away 
from  nominal  value  and  change  the  energy  state  of  the  aircraft.  This  energy  change  should 
be  lower  than  two  dimensional  maneuvers.  In  order  to  accomplish  this,  a  new  control 
law  governing  VW,  has  been  determined.  In  its  most  basic  sense,  the  control  law  has  two 
parts:  1)  energy  conservation  part  2)  energy  changing  part.  As  a  new  or  modified  way 
of  accomplishing  energy  excursion  minimizing  maneuvers,  allow  the  ew  term  in  Equation 
(5.23)  to  be  zero.  This  is  mathematically  correct,  since  evv  is  a  perturbation  quantity  and 
is  set  to  zero  to  obtain  the  energy  conserving  part  of  VW,.  Setting  the  derivative  of  a 
perturbation  to  zero  is  the  same  as  setting  the  perturbation  to  zero  because  a  perturbation 
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if  loo  at  time  equals  icro.  therefore,  the  wing  aircraft  specific  energy  perturbaticm,  ew, 
must  equal  sero.  In  equation  form, 


VWc  =  “iVir  [r - ”“1  W  -  (- — W.)  +  +  kg,iy  f  h,dt  (5.64) 

\  ITlkw  ^wJ  /  Jo 

This  is  now  implemented  to  check  the  results.  The  state  space  representation  of  the  system 
including  this  new  Vw,  is  listed  in  Appendix  D. 

5. 7  Dimensional  Eigenvalue  Analj/ais  For  New  Vwc 

(See  Appendix  A  for  ditnenaional  plant  models)  Like  Section  5.5,  the  objective  of 
this  section  is  to  compare  the  new  Laplace  transform  of  the  step  response  with  Busogany’s 
response  [2:page  6-5].  However,  in  this  analysis  energy  excursion  minimising  with  no  ew  in 
is  used.  In  order  to  do  this,  the  output  matrix,  Cg,  has  been  been  changed  to  reflect 
Busogany’s  output,  and  thus  examine  the  wing  velocity  response  (See  Section  5.5) 

The  dimensional  wing  velocity’s  step  response  is: 


yw(s)  =  fo  -  Ko  0  0  ol  [s/  -  r„  - 

L  ro  J  * 

- - - - ' 

Dimensional 

Ti-  /  ^  1  ,  ’’a  ,  .  **4  ,  ’’s  ,  *’6  . 

Vlv(s)  =  -  -I - ; - f-  — ; - 1 - ; - 1 - ; - 1 - ; - -| - ; - 

s  s+Pa  S  +  P3  «  +  P4  s  +  Ps  s-l-Ps  s-l-pr 


(5.65) 

(5.66) 


As  in  Section  5.5,  it  is  desired  to  perform  an  eigenvalue  analysis  to  determine  if  a 
fast  response  is  possible  for  Vw  This  is  accomplished  experimentally  by  keeping  the  ratio 
of  =  100  while  varying  their  respective  values.  This  is  shown  in  Figure  5.10. 

Figure  5.10  yields  some  interesting  observations.  In  order  to  track  a  step,  must 
be  >  .01.  kx,pv  must  also  be  <  1  so  that  p4  and  p5  do  not  go  unstable.  The  fastest 
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pole  it  p7.  It  remain*  at  -3  until  k,tpy  =  .1,  however,  its  residue  is  very  snudl  (»  0)  until 
ka,py  =  Yet,  its  residue  becomes  approximately  -.1  at  k^apy  =  1.  This  is  still  5  times 
less  than  the  residue  for  poles  4  and  5.  Its  contribution  would  be  negligible  and  would  make 
poles  4  and  5  slower  (i.e.  closer  to  the  imaginary  axis).  Therefore,  the  fastest  poles  with  the 
possibility  of  having  a  dominant  response  are  poles  4  and  5,  which  are  complex  coiyugates. 
They  are  the  fastest  poles  with  the  largest  residues.  Therefore,  picking  k^^py  =  .01  and 
kgMiv  -0001  results  in 


•  Vw  is  able  to  track  step  inputs 

•  Poles  p2,  p3,  and  p6  have  small  residues  they  do  not  contribute  much  to  the  output 
response 

•  Poles  p4  and  p5  are  complex  coiyug^t^  of  other  and  are  dominant 

•  Pole  p7  is  the  fastest  pole  but  its  residue  is  never  large  enough  to  have  an  impact  on 
the  output  response 

This  chapter  has  explored  two  ways  of  energy  minimising  maneuvers.  One  had  in 
the  formulation  for  VW<,i  while  the  other  disregarded  cw  Comparing  eigenvalue  analyses 
of  both  formiilations  will  show  the  differences. 

Transfer  function  and  Partial  FVaction  expansion  of  the  step  response  of  system  without 
ew- 


Vw{*) 

Vl{s) 

Vw{») 


(5.67) 


.1287(s  +  .5)(s  +  3.3486)(s  +  .0010441)(s  +  .01) 

(s  +  .5)(s  +  .24136  ±  i.28202)(s  +  .01041)(s  +  .0010437)(s  +  3.0058) 

1  -4.2185x10"^  .042551  -.52011  -  ;.45447  -.52011  +  7.45447 

-  + - .  .  +  ■  ' 


s  +  .0010437  «  +  .01041  s  +  .24136  +  j.28202  s  +  .24136  -  j.28202 


Dominant  Poles 


s  +  .5 


+ 


-1.9009  X  10-^ 
s  +  3.0058 


(5.68) 
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Examining  Equations  (5.62)  and  (5.67),  shows  Vw.  without  in  its  formulation 
should  yield  superior  performance  due  to  the  fact  that  the  system  without  ew  will  respond 
quicker  than  the  other  system  since  its  dominant  poles  are  fasta  (i.e.,  farther  to  the  left  in 
the  Real-Imaginary  Plane).  In  addition,  the  poles  in  Equation  (5.68)  are  faster  overaU  than 
those  of  Equation  (5.67),  and  Equation  (5.68)  has  only  one  pair  of  complex  poles  versus 
two.  The  extra  pair  of  complex  poles  in  Equation  (5.63)  are  very  slow  (i.e.  very  close  to 
the  imaginary  axis)  and  their  residues  are  larger  than  the  residues  of  the  corresponding 
real  poles  (p2  and  p7).  Since  the  denominator  has  3  more  poles  than  seros,  the  residues 
should  and  do  add  up  to  0  [5:page  117]. 

It  is  interesting  to  note  that  the  newer  formulation  for  Vw,  has  poles  close  to  Buso- 
gany’s  response.  Also,  leaving  ew  out  of  Vw^  allows  the  integral  and  proportional  gain  to 
be  increased. 


5.7.1  Nonlinear  Simidations  With  New  Vw^-  The  gains  for  the  heading  changes 
are  listed  in  Table  5.3,  and  the  responses  are  in  Figures  5.11  and  5.12.  The  gains  were 
determined  in  the  previous  section.  The  heading  changes  are  very  good.  The  responses 
are  almost  as  good  as  in  energy  conserving  maneuvers.  There  is  a  slightly  worse  response 
in  X  separation.  However,  this  response  is  much  better  than  previously  shown  when  ew 
was  included  in  the  feedback  control  law.  A  90°  heading  change  using  energy  conserv¬ 
ing/excursion  minimising  maneuvers  using  transposition  is  depicted  in  Figure  5.13.  The 
transposition  allows  lower  energy  usage  during  large  heading  changes. 


Table  5.3  PI  Controller  Values  For  Equation  (5.64)  and  Nonlinear  Simulations  Using  New 
Vw...  Heading  Changes  In  Figures  5.11  and  5.12 


Parameter 

Dimensional  Value 

^**PV 

.01 

.0001 

The  gains  for  the  velocity  change  are  given  in  Table  5.4.  The  25  velocity  command 
response  is  shown  in  Figure  5.14. 
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Thble  5.4  PI  CimtroUar  Values  For  Equation  (5.64)  and  Nonlinear  Simulation  Using  New 
Vw,:  Velocity  Change  In  Figure  5.14 


Parameter 

Dimensional  Value 

.02 

mam 

.002 

An  interesting  observation  can  be  discerned  concerning  the  velocity  change.  The 
proportional  and  integral  gains  cannot  be  chosen  such  that  the  ratio  between  proportional 
and  integral  gains  equals  100,  if  a  quick  response  is  desired.  The  integral  gain  had  to  be 
increased  in  order  that  the  recovery  from  the  altitude  change  would  occur  in  a  reascmable 
amount  of  time.  Otherwise,  recovery  is  on  the  order  of  600  -  800  seconds.  The  values 
in  Table  5.4  are  determined  heuristically.  Using  the  gains  determined  in  Section  5.7,  the 
gains  were  modified  in  a  trial  and  error  process  until  a  desirable  response  was  attained. 
Increasing  integral  gain  allows  the  altitude  to  respond  much  quicker. 

In  addition,  altitude  responses  are  investigated  in  this  section.  Like  velocity  changes, 
Busogany’s  formulation  could  not  handle  altitude  changes.  Busogany’s  control  in  energy 
conservation  consisted  of  controlling  altitude  through  z-separation  and  velocity  through 
altitude,  ff  the  lead  aircraft  changed  altitude,  there  was  no  way  for  the  formation  to  follow 
using  energy  conservation  because  there  was  no  controller  dependent  on  altitude  error. 
The  new  formulation  addresses  this  through  the  PI  controller  augmented  to  the  Vw,  iu 
Equation  (5.64).  This  controller  operates  on  altitude  error  and  thus  increases  the  speed 
of  the  wing  by  throttle  usage  when  there  is  an  altitude  differential.  This  is  seen  in  Figure 
5.15.  The  gains  used  for  the  PI  controller  in  Vw,  are  larger  than  those  in  Table  5.4.  They 
are  listed  in  Table  5.5. 

Table  5.5  PI  Controller  Values  For  Equation  (5.64)  and  Nonlinear  Simulation  Using  New 
Vw,:  Altitude  Change  In  Figure  5.15 


Parameter 

Dimensional  Value 

^**Pv 

.07 

mmm 

.007 
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Af  in  the  velocity  reepcmae,  the  int^rai  gain  had  to  be  increaaed  in  (xder  to  allow 
recovery  Cram  an  altitude  change  in  100  tecaadi. 

5.8  Conchuions 

This  ch^ter  has  introduced  the  concept  of  energy  excnrskm  tninimiMitig  maneuvers. 
It  explored  two  methods  of  minimising  energy  excursion  during  maneuvers.  Based  on  the 
results  presented,  the  wing  velocity  control  law  that  sets  ei^  =  0  results  in  better  responses 
(See  Equaticm  (5.64)).  Therefore,  it  is  the  preferred  method  of  achieving  energy  excurskm 
minimising  maneuvers. 
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Figure  5.13  Diamond  Formation  90°  Heading  and  Formation  Change  Utilising  Energy 
Consenring/Excursion  Minimising  [2;Figure  6.6] 
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VI.  Comparison  of  Two  Dimensional  and  Three  Dimensional  Maneuvers 

Two  Dimensional  Maneuvers  and  Three  Dimensional  Maneuvers  are  explored  in 
Chapters  IV  and  V,  respectively.  The  objective  of  this  chapter  is  to  simulate  and  compare 
the  following  methods: 

1.  Two  dimensional  maneuvers  (i.e.,  all  maneuvers  occur  in  the  same  plane  of  flight) 

2.  Energy  conserving  maneuvers  (i.e.,  all  maneuvers  occtir  with  constant  aircraft  specific 
energy,  in  and  out  of  the  same  plane  of  flight) 

3.  Energy  excursion  minimising  maneuvers  (i.e,  all  maneuvers  occur  without  constant 
aircraft  specific  energy,  in  and  out  of  the  same  plane  of  fiight) 

The  results  of  this  comparison  are  evaluated  and  a  recommendation  for  further  study 
is  stated. 

6. 1  Nonlinear  Simulations 

This  section  compares  the  responses  for  the  nonlinear  simulations  of  the  three  meth¬ 
ods  outlined  in  the  introduction  of  this  chapter.  The  gains  used  are  listed  in  Table  6.1. 
Once  again,  these  are  the  gains  used  by  Buzogany  in  Chapters  5  and  6  of  his  thesis.  This 
is  done  to  show  the  eifec*  of  including  the  PI  controller  into  the  of  energy  minimising 
maneuvers.  In  all  the  plots,  solid  line  indicates  energy  conserving  maneuvers,  duhed  line 
indicates  energy  excursion  minimizing  maneuvers,  and  dotted  line  indicates  two  dimen¬ 
sional  maneuvers. 

6.S  Discussion 

Figures  6.1  -  6.4  yield  the  following  observations: 

1.  The  two  dimensional  maneuvers  provide  the  best  response  in  x-separation  overall. 
The  deviations  away  from  nominal  are  smaller  than  the  energy  conserving  and  energy 
excursion  minimizing  maneuvers 
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Figure  6.3  C-130A:  Diamond  Formation,  25  ^  Velocity  Change.  Nonlinear  Simu¬ 

lation:  Energy  Conserving  Implonentation  (solid  line),  Energy  Minimizing 
Implementation  (dashed  line),  and  Two  Dimensional  Implementation  (dotted 
line) 
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Figure  6.4  C-130A:  Nonlinear  Simulation,  Response  To  A  100  ft  Step  Altitude  Change: 

Energy  Conserving  Implementation  (solid  line),  Energy  Mifiimiging  Imple¬ 
mentation  (dashed  line),  and  Two  Dimensional  Implementation  (dotted  line) 
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Tbble  6.1  Gaini  Used  In  Nonlinear  Simulations  To  Compare  All  Three  Methods  For 
Heading,  Velocity,  and  Altitude  Changes  In  Figures  6.1  -  6.4 


Parameter 

Dimensional  Value 

Two  Dimensional  Maneuvers 

Energy  Conserving 

Energy  Minimising 

10*  a 

90*  a 

Valodlf  A 

AltUB4c  A 

AU  a> 

B««4img  A 

V^cilf  A 

AllltsAa  A 

^MZPV 

.01 

.02 

.07 

.0001 

.002 

.007 

k. 

1 

1 

1 

1 

1 

1 

1 

1 

ky 

0 

3 

0 

0 

0 

0 

0 

0 

kxp 

.75 

1 

.75 

.75 

.75 

.75 

.75 

.75 

k,i 

.0075 

.1 

.0075 

.0075 

.0075 

.0075 

.0075 

.0075 

fh 

1 

1 

1 

1 

1 

1 

1 

1 

k* 

■Qm 

4 

0 

0 

0 

0 

mom 

mjnn 

1 

1 

1 

1 

1 

1 

1 

1 

kpi 

.01 

.1 

.01 

.01 

.01 

.01 

.01 

.01 

k,p 

1 

.5 

2.  There  are  no  wing  aircraft  altitude  changes  using  two  dimenaional  maneuvers  for 
heading  and  velocity  changes 

3.  Two  dimensional  maneuvers  cause  the  largest  amount  of  wing  specific  energy 
excursions  =  throttle  usage 

4.  Energy  conserving  maneuvers  allow  heading  changes  that  have  a  worse  response 
in  x-se^  ^A'ation  (i.e.,  there  are  larger  deviations  away  from  nominal) 

5.  Energy  conserving  maneuvers  use  the  least  amount  of  energy  during  heading 
changes.  The  specific  energy  of  the  wing  aircraft  is  not  exactly  constant.  Due  to 
rate  limiters  in  the  nonlinear  simulation,  the  specific  energy  of  the  wing  aircraft  do 
not  remain  at  their  nominal  value.  However,  the  excursions  away  from  nominal  are 
small  and  negative 


6.  Energy  conserving  maneuvers  have  a  worse  x-separation  response.  They  also 
have  altitude  drops  associated  with  their  accomplishment.  In  addition,  the  velocity 
perturbs  farther  away  from  nominal  than  in  two  dimensional  maneuvers 
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7.  Energy  excursion  mininriiaing  maneuvers  cause  the  second  largest  amount  of 
energy  excursions  during  maneuvers.  This  is  expected.  It  should  be  between  energy 
conserving  maneuvers  and  two  dimensional  maneuvers. 

8.  Energy  excursion  minimising  maneuvers  also  yield  x-separation  performance 
between  two  dimensional  and  energy  conserving  maneuvers.  Additionally,  the  veloc¬ 
ity  response  shows  much  larger  excursions  away  firom  nominal  during  maneuvers 

9.  Energy  excursion  minimising  maneuvers  largest  advantskge  over  energy 
conserving  maneuvers  is  its  recovery  firom  velocity  and  altitude  changes. 
It  also  has  a  smaller  altitude  change  than  energy  conserving  maneuvers 

Observations  1-9  show  that  there  are  distinct  differences  between  the  three  methods. 
Two  dimensional  maneuvers  are  good  if  it  is  unwise  for  wing  aircraft  in  the  formation  to 
track  the  lead  aircraft  without  altitude  deviations  and  large  x-separation  deviations.  How¬ 
ever,  these  maneuvers  cause  significant  energy  excursions  more  so  than  other  maneuvers, 
resulting  in  more  fuel  use  by  the  wing  aircraft.  Energy  conserving  maneuvers  are  good  if 
only  interested  in  heading  changes.  Wing  aircraft  would  use  less  fuel  and  would  have  the 
same  energy  consumption  as  the  lead  aircraft.  Energy  excursion  minimising  maneuvers 
provide  the  best  overall  capabilities  and  energy  usage  if  altitude  changes  do  not  cause  any 
problem  within  the  formation. 

There  are  times  when  each  of  these  methods  is  appropriate  for  a  formation.  Two 
dimensional  maneuvers  use  the  most  energy.  However,  the  characteristics  of  no  altitude 
changes  and  smaller  deviations  in  x-separation  would  benefit  a  low  altitude  formation 
penetrating  hostile  territory. 

Both  energy  conserving  maneuvers  and  energy  excursion  minimising  maneuvers  pro¬ 
vide  good  characteristics  for  a  formation  in  friendly  skies  at  cruising  altitude.  Altitude 
changes  of  a  few  hundred  feet  are  of  no  consequence  at  20,000  feet  above  ground  level 
(AGL).  Separation  requirements  can  also  be  relaxed  since  detection  by  enemy  forces  is 
not  a  threat.  The  question  is  what  kind  of  performance  is  desired.  Heading  changes  are 
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definitely  better  when  done  with  energy  conserving  maneuvers.  The  aircraft  energy  is 
preserved  and  fuel  use  is  kept  consistent  with  the  lead  aircraft.  However,  velocity  changes 
cannot  be  performed  using  constant  energy  numeuvm.  Energy  excursicm  minimiging  ma¬ 
neuvers  are  capable  of  performing  velocity  changes  with  less  energy  consumption  than  two 
dimensionid  maneuvers.  Their  only  disadvantage  is  they  are  not  as  efficient  as  energy 
conserving  maneuvers  during  heading  changes.  If  one  or  the  other  must  be  chosen,  energy 
excursion  minitnivmg  nuneuvers  definitely  provide  the  most  capabilities  with  lower  use  of 
energy  than  two  dimensional  maneuvers. 

A  logical  extension  of  these  observations  can  be  determined  by  examining  the  ”fliight 
phases”  of  aircraft.  Nelson  divides  the  flight  phases  of  an  aircraft  into  the  foUowing: 

1.  Nonterminal  Flight  Phase 

(a)  Category  A  :  "Nonterminal  flight  phases  that  require  rapid  maneuvering,  pre¬ 
cision  tracking,  or  precise  flight-path  control.  Included  in  the  category  are 
air-to-mr  combat  ground  attack,  weapon  delivery /launch,  aerial  recovery,  recon¬ 
naissance,  in-flight  refueling  (receiver),  terrain-following,  antisubmarine  search, 
and  close-formation  flying” 

(b)  Category  B  :  "Nonterminal  flight  phases  that  are  normally  accomplished  using 
gradual  maneuvers  and  without  precision  tracking,  although  accurate  flight-path 
control  may  be  required.  Included  in  the  category  are  climb,  cruise,  loiter,  in¬ 
flight  refueling  (tanker),  descent,  emergency  descent,  emergency  deceleration, 
and  aerial  delivery” 

2.  Terminal  Flight  Phase 

(a)  Category  C  :  "Terminal  flight  phases  are  normally  accomplished  using  gradual 
maneuvers  and  usually  require  accurate  flight-path  control.  Included  in  this  cat¬ 
egory  are  takeoflf,  catapult  takeoff,  approach,  wave-off/go-aroimd,  and  landing” 

[8] 
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These  ddmiticms  of  flight  phase  concern  individual  aircraft.  However,  they  can  be  extended 
to  a  formation  of  aircraft  also.  Two  dimensional  maneuvers  would  be  ideal  for  Category 
A  flying.  This  application  could  very  easily  be  used  for  reconnaissance,  terrain-following, 
and  close-formation  flying.  Category  B  could  be  handled  by  energy  excursion  miniTniaing 
maneuvers  especially  cruising,  loitering,  climb,  or  descent.  Obviously,  the  appropriate  form 
of  control  would  need  to  be  used  for  the  corresponding  flight  phase.  Category  C  flight  phase 
will  probably  never  be  used  with  this  kind  of  control  system. 

6.S  Conchuion 

All  three  maneuvering  methods  compared  in  this  chapter  are  useful  concepts  of  auto¬ 
matic  formation  flight  control,  however,  each  has  its  advantage  and  disadvantage.  Different 
flight  phases  experienced  by  the  formation  would  determine  which  concept  to  use. 


6-9 


VII.  Analyaia  and  Conchuiona 
7. 1  Objectivea  of  Reaearck  Met 

The  objective  of  this  research  is  to  continue  the  previous  work  acannplished  by  Rohs, 
Dargan,  and  Busogany,  by  «^h«nring  the  Automatic  Formation  Flight  Ccmtrol  Syston  with 
the  inclusion  of  a  linear  mixer,  which  blends  velocity  and  heading  rate  with  the  longitudinal 
and  lateral  formation  spacing,  and  investigating  energy  excursion  minitniKitig  maneuvers. 
These  objectives  have  been  accomplished. 

7.1.1  Horizontal  Formation  Fligki  Control  -  Analyaia.  The  inclusion  of  a  linear 
mixer  into  the  Automated  Flight  Control  System  is  accomplished  in  Ch^ter  IV  for  both 
the  X  and  y  channels.  Two  dimensional  maneuvers  are  considered.  The  difference  between 
the  previous  effort  accomplished  by  Busogany  [2]  and  the  new  effort  is  the  inclusion  of  the 
linear  mixer  gains  and  lead  aircraft  states  in  the  analysis.  Busogany’s  analysis  fed  back 
only  separation  errors  while  the  new  formulation  feeds  back  separation  errors,  heading  error 
(V*«,  y-channel):  and  velocity  error  (V«,  x-channel).  Previously,  the  lead  aircraft  states  were 
considered  to  be  disturbances  into  the  system.  Now,  the  disturbances  into  the  system  are 
the  lead  aircraft  commands.  The  lead  aircraft  states  augment  the  plant’s  state  vector.  Be¬ 
cause  the  six  state  model  is  used  with  first-order  aircraft  chaimel  models,  the  characteristic 
polynomials  for  each  channel  are  of  fourth-order.  Performing  a  Routhian  stability  analysis 
on  each  channel  results  in  conditions  for  system  stability.  A  gain  parameter  stability  enve¬ 
lope  in  three  dimensional  space  is  determined  by  performing  a  linear  transformation  on  the 
gain  parameters  of  interest.  The  gains  of  interest  for  each  channel  are:  k^,  and  k^  for 
the  y-channel  and  k,,,  and  k,  for  the  c-channel  (recall:  k,  =  kg  =  1).  There  are  some 
quadratic  terms  in  each  channel  that  make  plotting  the  envelope  difficult.  By  analysing 
the  stability  equations,  a  nonlinear  transformation  of  variables  is  determined  between  the 
e  and  Gy  =  (*,/,*^)  €  n\  and  Kx  =  6  71^  and 

Gx  =  (nt,  n,  kgi)  parameter  vectors.  These  nonlinear  transformations  map  the  respective 
gains  into  a  space  in  which  the  stability  equations  become  more  tractable  and  thus  easier 
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to  plot.  With  these  stability  oivelopes,  any  choice  of  gains  can  be  transformed  into  the 
new  space  to  graphically  gauge  the  stability  of  the  system. 

Simulations  with  the  mixer  in  place  showed  better  response  characteristics  in  ve¬ 
locity  and  separation  error.  Thus,  including  mixers  and  lead  aircraft  states  into  the  new 
Formation  Flight  Control  model  has  successfully  been  completed  and  validated. 

1.1. i  Three  Dimetuional  Energy  Exatnion  Minimizing  Maneuvers.  Previously, 
energy  conserving  maneuvers  were  accomplished  by  keeping  the  specific  energy  of  the  wing 
aircraft  constant.  Thus,  any  increase  in  lead  aircraft  velocity  results  in  a  corresponding 
decrease  in  wing  aircraft  altitude,  and  vice-versa.  Hence,  energy  conserving  maneuvers  are 
fine  during  formation  heading  changes,  however,  velocity  or  altitude  changes  result  in  loss 
of  formation  integrity.  This  research  effort  allows  the  specific  energy  of  the  wing  aircraft 
to  deviate  from  its  nominal  value  by  including  a  PI  controller  for  Vw,-  Thus,  a  PI  con¬ 
troller  commands  the  throttle  (Vw.)  which  causes  excursions  in  specific  energy.  Therefore, 
formation  energy  changes  can  be  made  and  formation  integrity  will  be  maintained. 

7.1.S  Comparison  Of  Formulations.  The  three  formulations  for  Automatic  For¬ 
mation  investigated  in  this  thesis  are: 

•  Two  dimensional  formation  maneuvers 

•  Energy  conserving  maneuvers  (i.e.,  specific  energy  is  kept  constant) 

•  Energy  excursion  minimising  maneuvers  (i.e.,  specific  energy  is  allowed  to  make 

excursions  away  from  its  nominal  value) 

As  expected,  the  energy  conserving  maneuvers  provided  heading  changes  with  mir.imal  ex¬ 
cursions  in  specific  energy  due  to  aircraft  rate  limiters.  However,  they  do  not  provide  good 
response  characteristics  for  velocity  changes  or  altitude  changes.  Both  two  dimensional 
and  energy  excursion  minimising  maneuvers  provide  formation  energy  change  caj^^abilities 
in  addition  to  heading  change  capability.  The  two  dimensional  maneuvers  provide  better 
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performance  overall  in  z  separation,  y  separation,  wing  velocity,  and  wing  altitude.  How¬ 
ever,  this  is  done  at  the  cost  of  requiring  more  wing  aircraft  specific  raergy  excursions  than 
the  energy  minimising  maneuvers,  which  change  the  altitude  in  order  to  compensate  for 
velocity  changes. 

These  formulations  for  Automatic  Flight  Control  provide  different  degrees  of  per¬ 
formance  and  energy  use.  There  is  a  tradeoff  between  capability  and  energy  use.  This 
is  a  useful  finding  because  it  forms  the  foundation  for  a  control  system  that  provides  de¬ 
sirable  characteristics  in  Flight  Phases  A  -  C  [8].  In  cruise,  a  blend  of  energy  conserving 
and  energy  excursion  minimising  would  provide  sufficient  formation  capabilities.  However, 
penetrating  enemy  airspace  would  obviate  the  need  for  two  dimensional  maneuvers,  where 
tighter  formations  and  lower  altitudes  could  be  flown. 
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-i-  =  best  response 
o  =  good  response 
—  =  worst  response 
n/c  =  not  capable 

1.2  Conclusions  and  Lessons  Learned 

The  following  conclusions  are  drawn  from  this  research  effort: 

•  In  2iccordance  with  research  accomplished  by  Busogany  in  December,  1992,  analysis 
proves  the  necessary  and  sufficient  condition  of  integral  control  action  to  achieve  zero 
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steady-state  error  for  a  formation  of  aircraft.  Additionally,  a  parameter  stability 
envelope  can  be  plotted  to  show  regions  of  guaranteed  stability  for  the  gains  o{ 
interest. 

•  Using  the  linearised,  state  space  MIMO  plant  model,  PI  controller  gains  can  be 
determined  through  pole  placement  techniques.  These  gains  provide  good  response 
characteristics  in  the  linear  case,  however,  they  need  to  be  adjusted  using  a  heuristic 
process,  due  to  the  rate  limiters  in  the  nonlinear  aircraft  channel  models. 

•  Adding  a  linear  mixer,  which  blends  velocity  and  heading  rate  with  formation  spac¬ 
ing,  into  the  analysis  yields  a  better  overall  response. 

•  Analysis  of  energy  excursion  minimising  maneuvers  also  proved  the  necessary  and 
sufficient  condition  of  integral  control  action.  Due  to  time  constraints,  stability 
envelopes  were  not  researched.  However,  the  linearised  MIMO  plant  model  allowed 
the  use  of  pole  placement  techniques  to  provide  values  of  the  PI  controller  parameters. 
Once  again,  a  heuristic/experimentation  process  is  required  to  determine  the  best 
gains  for  the  full  nonlinear  simulation. 

•  Comparison  of  the  three  concepts  for  Automatic  Formation  Flight  Control  show  the 
relative  capabilities  of  each  method.  These  capabilities  can  be  used  to  synthesise  a 
flexible  control  system  optimised  for  different  phases  of  flight. 

7.3  Recommendations  For  Further  Study 

As  stated  in  Chapter  I,  this  thesis  is  a  stepping  stone.  Some  recommendations  for 
further  study  are; 

1.  Include  sensor  models  into  the  simulation.  Additionally,  sensor  noise  should  be  in¬ 
corporated  to  evaluate  its  effect  on  system  stability  [2;page  8-5]. 

2.  Investigate  incorporating  air-refueling  capabilities  into  Automatic  Formation  Flight 
Control 
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3.  Determine  a  method  to  quantify  the  effect  of  nonlinear  rate  limiters  oa  the  control 
system  response  in  order  to  remove  the  heuristic  process  of  selecting  “good”  controller 
and  mixer  gains 

4.  Include  second  ord»  models  developed  by  Busogany  and  perform  analytic  analyses 
to  determine  the  pertinent  parametric  stability  envelopes 

5.  Time  delays  representing  update  data  rates/sampling  rates  should  be  incorporated 
into  the  model  and  simulation  in  order  to  quantify  the  effects  on  system  stability 
[4:page  6-4] 

6.  Investigate  energy  conserving  and  excursion  Tnifiifni«ing  maneuvers  with  disaimilar 
aircraft 

7.  Energy  excursion  miniTnising  maneuvers  should  be  utilized  with  multiple  wing  air¬ 
craft  to  investigate  its  effect  on  formation  integrity 

8.  A  robustness  investigation  applying  the  Routh  Stability  Criterion  for  both  two  di¬ 
mensional  and  energy  minimizing  maneuvers  should  be  explored 

9.  Optimal  control  techniques  should  be  applied  to  energy  minimization  maneuvers  in 
order  that  the  optimal  energy  use  during  maneuvers  is  determined 

This  concept  can  also  be  applied  to  other  vehicles  that  move  in  large  groups.  Some 
future  applications  are: 

1.  Extend  this  research  to  helicopters.  Rotorblades  make  flying  in  tight  formation 
dangerous 

2.  Explore  formations  composed  of  fighter  aircraft,  to  explore  high  performance  consid¬ 
erations 

3.  Explore  formations  composed  of  hovercraft,  which  are  used  for  beach  assaults 
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7.^  Summary 

As  the  research  performed  by  Biuogany  and  Dargan  showed  [2]  [4],  a  PI  controUer  can 
successfully  maintain  a  multiple  aircraft  formation  through  different  maneuvers.  Analysis 
reaflirmed  the  necessary  and  sufficient  condition  of  integral  control  action  for  maintaining 
formation  integrity  during  maneuvers.  Parameter  stability  envelopes  can  be  determined  in 
the  two  dimensional  case.  Two  dimensional  and  energy  excursion  minimising  maneuvers 
provide  desirable  formation  response  characteristics  for  the  appropriate  flight  phase. 
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Appendix  A.  Dimensional  Plant  Matrices 


Through  the  course  of  this  thesis,  analysis  is  performed  on  the  linearised  state  space 
MIMO  model  and  the  dimensional  eigenvalues.  Chapters  IV  and  V  interchange  between 
dimensional  and  nondimensional  discussions.  In  order  to  avoid  confusion,  the  MIMO  state 
space  models  in  these  chapters  are  presented  in  their  dimensional  form 


A.l  Chapter  IV  DimensiontU  Models 
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A. 2  Chapter  V  Dimenaional  Model 
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Appendix  B.  Stability  Equations  Prom  Chapter  V 


The  following  stability  equations  are  from  Section  5.1.3.  They  are  derived  from  a 
Routhian  array  in  the  same  section.  Due  to  the  complexity  of  the  algebra,  Mathematica 
was  used  to  derive  the  larger  inequalities  in  this  appendix.  For  more  information,  see 
Section  5.1.3. 
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Appendix  C.  Y  and  X  Channel  Stahmty  Conatraints  PloUed  Separately 

The  stability  constraints  of  Equations  (4.84)  -  (4.87)  for  the  y-channel  and  Equatiosu 
(4.88)  •  (4.91)  for  the  c-channel  are  plotted  individually  and  then  together  on  the  same  plot 
in  this  appendix.  Figures  4.3  and  4.5  show  the  stability  envel<^>es  of  the  y  and  c-channel, 
respectively.  This  iq>pendix  shows  all  of  the  intermediate  plots  that  were  used  to  derive 
the  final  stability  envelopes. 

•  Figures  C.l  -  C.5  and  Figures  C.7  -  C.ll  show  the  intermediate  plots  of  the  pertinent 
stability  constraints.  The  corresponding  equation  number  is  listed  in  the  caption  of 
each  figure 

•  Figures  C.6  and  C.12  show  the  final  stability  envelopes  of  the  y  and  c-channels,  as  in 
Chapter  IV.  They  are  included  here  to  show  the  complete  process,  from  beginning 
to  end 
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Figure  C.l  Plot  of  the  Nondimensional  Y-Channel  Stability  Constraint  for  Rx>w  s’  of  the 
Routhian  Array  (Equation  (4.84)) 
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Figure  C.2  Plot  of  the  Nondimensional  Y-Channel  Stability  Constraint  for  Row  s’  of  the 
Routhian  Array  (Equation  (4.85)) 
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Figure  C.3  Plot  of  the  Nondimensional  Y-Channel  Stability  Constraint  for  Row  of  the 
Routhian  Array  (Elquation  (4.86)) 
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Figure  C.6  Nondimensional  Y-Channel  Envelope;  Stability  Between  the  Surfaces 
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Figure  C.ll  All  of  the  Nondimensional  X- Channel  Constraints 


The  Stability  Envelope  for  die  X-Channel  from  Aoodier  Angle 
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Figure  C.12  Nondimeiuional  Y-Channel  Envelope;  Stability  Between  the  Surfaces 
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Appendix  D.  Energy  Minimizing  Maneuvers  WMiout  ew  in  VW< 


In  Chapter  V,  it  was  determined  Vw,  formulated  without  ew  provided  better  tran¬ 
sient  time  responses  and  energy  usage  than  the  system  including  ew-  This  is  attributed  to 
the  assumption  that  there  are  two  different  parts  of  Vw/-  an  energy  conserving  part  and 
an  energy  using  part.  Together,  they  yield  an  energy  minimising  control  law.  The  energy 
conserving  portion  of  the  control  law  is  derived  with  the  assumption  that  the  specific  en¬ 
ergy  perturbation,  ew,  is  sero.  However,  when  the  PI  controller  is  added  to  this  energy 
conserving  control  law,  the  assumption  that  ew  remains  at  sero  is  violated  since  the  PI 
controller  is  directly  driving  the  ew  perturbation.  Thus,  the  resulting  model  is: 
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and, 

=  augmented  xz-channel  state  vector 
=  augmented  xz-channel  plant  matrix 
=  augmented  xz-channel  disturbance  input  matrix 
=  augmented  xz-channel  disturbance  input  vector 
=  augmented  cz-channel  output  vector 
=  augmented  xz-channel  output  matrix 

The  only  difference  between  this  model  and  the  one  in  Chapter  V  is  the  coefficients 
Dkvi  ^twi  ^hti  find  These  are  the  terms  which  contain  those  expressions  related 

to  the  derivative  ofew 

D.  1  Static  Stability  Analysis 

A  static  stability  analysis  has  been  performed  upon  the  system  with  no  ew  included 
in  the  Vw,-  This  analysis  is  the  same  as  in  Section  5.1.2.  The  determined  result  is  the 
same  answer  that  the  previous  analysis  yielded.  See  Equation  (5.39). 

D.2  Dynamic  Stability  Analysis 
The  characteristic  equation  is; 

det(5/  -  )  =  z*  -1-  ps*  qs*  -|-  rz’  -|-  tz*  -f-  us  -|-  v  (D.7) 

where  the  coefficients  are, 
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The  Routhian  array  resulting  from  this  characteristic  equation  is  shown  in  Table  D.l. 
Table  D.l  XZ- Channel  Routhian  Array 
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The  resulting  stability  equations  resulting  from  applying  the  Routh  Stability  Crite¬ 
rion  to  Table  D.l  are: 
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row  j* 
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rows  are  too  complex  and  large  to  include 

row 
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In  this  case,  where  ew  is  absent  firom  Vw,,  integral  control  action  is  required  for  both 
zero  steady-state  error  and  transient  or  dynamic  stability. 
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The  research  contained  in  this  thesis  explores  the  concepts  of  Automated  Formation  Flight  Control  documented 
in  three  previous  AFIT  theses.  The  generic  formation  analyzed  consists  of  a  Leader  and  Wingman,  with  the 
Wingman  referencing  its  maneuvers  off  of  Leader  maneuvers.  Specifically,  planar  formation  flight  control  con¬ 
cerning  only  heading  and  velocity  changes  is  considered.  Next,  the  vertical  separation  constraint  is  relaxed  to 
allow  wing  maneuvers  outside  of  the  flight  plane  of  the  lead  in  order  to  minimize  the  energy  expended  by  the 
wing  in  a  maneuver.  Analysis  of  the  two  forms  of  formation  flight  control  investigated  in  this  thesis  reveals  the  > 
close  relationship  between  formation  geometry,  aircraft  time  constants,  controller  gains,  formation  performance,  ' 
and  control  system  stability.  Integral  control  action  is  determined  to  be  necessary  for  formation  flight  control,  i 
Nonlinear  simulations  are  accomplished  on  a  digital  computer  to  validate  the  analysis  of  the  automated  formation  i 
flight  control  system.  Compuisons  are  made  between  the  two  forms  of  formation  flight  control  considered,  and 
a  third,  energy  conserving  maneuvers,  in  order  to  determine  which  is  best  for  each  phase  of  flight. 
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